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Abstract 

Context. The low metallicity interstellar medium (ISM) is profoundly different from that of normal systems, being clumpy with low 
dust abundance and little CO-traced molecular gas. Yet many dwarf galaxies in the nearby universe are actively forming stars. As the 
complex ISM phases are spatially mixed with each other, detailed modeling is needed to understand the gas emission and subsequent 
composition and structure of the ISM. 

Aims. Our goal is to describe the multi-phase ISM of the infrared bright low-metallicity galaxy Haro 11, dissecting the photoionised 
and photodissociated gas components. 

Methods. We present observations of the mid-infrared and far-infrared fine-structure cooling lines obtained with the SpitzerflRS and 
HerschelfPACS spectrometers. We use the spectral synthesis code Cloudy to methodically model the ionised and neutral gas from 
which these lines originate. 

Results. We find that the mid- and far-infrared lines account for ~1% of the total infrared luminosity Ltir, acting as major coolants of 
the gas. Haro 1 1 is undergoing a phase of intense star formation, as traced by the brightest line, [O in] 88 //m, with L [0 hij/Ltir ~ 0.3%, 
and high ratios of [Ne m]/[Ne n] and [S iv]/[S m]. Due to their different origins, the observed lines require a multi-phase modeling 
comprising: a compact H n region, dense fragmented photodissociation regions (PDRs), a diffuse extended low-ionisation/neutral gas 
which has a volume filling factor of at least 90%, and porous warm dust in proximity to the stellar source. For a more realistic picture 
of the ISM of Haro 1 1 we would need to model the clumpy source and gas structures. We combine these 4 model components to 
explain the emission of 17 spectral lines, investigate the global energy balance of the galaxy through its spectral energy distribution, 
and establish a phase mass inventory. While the ionic emission lines of Haro 1 1 essentially originate from the dense H n region 
component, a diffuse low-ionisation gas is needed to explain the [Ne n], [N n], and [C n] line intensities. The [O in] 88 /im line 
intensity is not fully reproduced by our model, hinting towards the possible presence of yet another low-density high-ionisation 
medium. The [O i] emission is consistent with a dense PDR of low covering factor, and we find no evidence for an X-ray dominated 
component. The PDR component accounts for only 10% of the [C n] emission. Magnetic fields, known to be strong in star-forming 
regions, may dominate the pressure in the PDR. For example, for field strengths of the order of 100 yuG, up to 50% of the [C n] 
emission may come from the PDR. 

Key words. galaxies:ISM - galaxies:individual (Haro 11) - ISM: line and bands - ISM: structure - techniques: spectroscopic - 
radiative transfer 



1. Introduction 

Star formation in the most pristine environments of the early 
universe is poorly understood. The closest analogs to chemi- 
cally unevolved systems are the low metallicity dwarf galax- 
ies of our local universe. But even those well-studied nearby 
galaxies continue to plague us with fundamental questions on 



* Herschel is an ESA space observatory with science instruments 
provided by European-led Principal Investigator consortia and with im- 
portant participation from NASA. 



what triggers and regulates star formation in metal-poor gas- 
rich conditions. Certainly, from an observational point of view, 
at all wavelengths, local universe low-metallicity dwarf galax- 
ies show dramatic differ ences comp ared to their more metal-rich 
counterparts dKunth & O stlin 2000, for a review), in the mid- 
infrared (MIR), far-infrared (FIR), and submillimeter (submm) 
wavelength regimes. These include the dearth of Poly cyclic 
Aromatic Hydrocarbons (PAHs), prominent hot MIR-emitting 
dust, submm excess, and enhanced [C n] 157 jum/CO(l-0) ratios 
(e.g. lMaddenll2008HMadden et aUl2012al) . 
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What physical properties, local as well as global, within the 
dwarf galaxies are we witnessing with these observational sig- 
natures? How do we turn these signatures into a realistic view 
of star formation under early universe conditions? How do basic 
local parameters, such as radiation field, density, compactness, 
filling factors, metallicity, etc. figure in the picture we have of 
star-forming dwarf galaxies? Much of the ambiguity surround- 
ing these questions is due to the lack of understanding the pre- 
cise role of critical diagnostics and of spatial resolution in the 
most nearby objects, especially in the FIR/submm, resulting in 
a mixture of different interstellar medium (ISM) phases in the 
integrated view of galaxies. 

Low metallicity star-f orming dwarf galaxies host young 
massive stellar clusters (e.g. Turner et al.l2000t[Beck et al. 2000; 



iJohnson & Kobulnickvl l2003) which produce ultraviolet (UV) 
photons. The low abundance and prevailing hard radiation field 
can explain most of the differences observed in dwarf galaxies 
in destroying large grains and PAHs, in favor of smaller dust 
grains. The UV photons also photodissociate molecules such 
as CO, which are unable to adequately self-shiel ding, thus re- 
ducing the CO abundance in metal-poor gas (e.g. IWolfire et al.l 
l2010tlNaravanan et alJl2012h . Thus UV photons play an im- 
portant role in the heating and the chemistry of the ISM, con- 
trolling the structure of H n regions and photodissociation re- 
gions (PDRs). While the heating of the gas is mainly due 
to photoionisation and gas-grain collisions in the H n regions 
and to the photoelectric effect and cosmic rays in the neu- 
tral phase, the cooling is normally dominated by radiative de- 
excitation of fine-structure lines, readily observed in the MIR 
and FIR/submm. The behavior of these tracers and their rela- 
tionship to the star formation properties, have been extensively 
studied to access physical parameters in terms of the gas tem- 
perature, density, and the radiation fie ld dMalhotra et al.l [200 1 ; 
iHunter et al.ll200ltlBrauher et al.ll2008l e.g.). These sudies make 
use of p hotoionisation and PDR models at the scale of individual 
cloud s dKaufman et al Ill999t IWolfire et al.lll990h iFerland et al 



1998: Abel et al.ll2005tlRollig et al.ll2006l: iTielens & Hollen bach 



19851) . 

The [C n] 157 fim line is an important FIR diagnostic of- 
ten put forward as a star formation trac er dStacev et al.1 120101: 
iBoselli et aDl2002t Ide Looze etai]l201ll) . It has been observed 
in man y different object ty pes, including dense star-forming re- 
gions (|Stacev et ak| J 199 lh . diffuse ionised and atomic regions 
(Madden et al. 1993), and high-redshift galaxies (Maiolino et al. 
2009: IStacev et al.ll2(")IO). Since the [C n] line is w ell-correlated 
with the 12 CO(1-0) line a nd the FIR luminosity (e.g. lStacev et al.1 
ll99lUPierini et al.l2 003). it is generally considered that the [C n] 
line traces massive star formation within the PDR paradigm. 
[C n] is normally the brightest FIR cooling line, being much 
brighter than CO, and carrying about 0.1 to ~0.5% of the FIR 
luminosity in metal-rich star-forming galaxies. In low metallic- 
ity dwarf galaxies, where CO( l-O) is not easily detected be- 
cause of photodissociation (e.g. lSchruba et alJ2012tlLerov et alJ 
1201 lh . [C n] can be as high as 1% of the L FIR and often 5 to 10 
times brighter relative to CO than in dustier starburst galaxies 
(e.g. IPoglitsch et al 1 1 1 9951: iMadden et all 1 19951: lMaddenll2000t 
IHunter et al.ll2001l: ICormier et al.ll2010t) The high [C n]/CO(l- 
0) ratio in low-metallicity systems may be alerting us to the 
fact that the morphology of their PDRs and H n regions and 
the relative filling factors of the various phases of dwa rf galax- 
ies differ from those of more metal-rich obj ects (e.g. iMaddenl 
2000; Kau fman et al.ll2006l:lR"oUig et alJ2006l) . Even in the near- 
est neighboring dwarf galaxies, the picture of star formation - 
and what controls this star formation - is enigmatic. 



However, using [C n] alone can lead to ambiguous interpre- 
tation, particularly on the origin of the C + emission, which can 
arise from the ionised as well as neutral phases of the ISM. To 
employ the [C n] line as a valuable star formation tracer, it is 
necessary to model the different phases on galaxy-wide scales. 
The challenge rests on the unavoidable fact that on the scale of 
galaxies, aside from the closest dwarf galaxies such as the LMC 
(50 kpc, 1/2 Z ), SMC (60 kpc, 1/5 Z ), NGC6822 (500 kpc, 
1/4 Z ), or IC 10 (800 kpc, 1/3 Z ), the phases of the ISM are 
mixed together in single telescope beams. When using a lim- 
ited number of diagnostic lines, only a single galactic compo- 
nent can be modeled and average conditions derived, which is 
otherwise an ensemble of individual PDRs, ionised regions, etc. 
This requires modeling a comprehensive suit of tracers covering 
a range of critical densities and excitation potentials to charac- 
terise the dense H n regions, the diffuse H n regions, PDRs, etc. 
self-consistently. This approach allow us to grapple with degen- 
eracies and ambiguities in the interpretation of these tracers. In 
this way, valuable observational diagnostics can be turned into a 
more accurate description of the multi-phase configuration of a 
full galaxy. 

The Herschel Space Observatory (Pilb ratt et al.l 12010) has 
brought the opportunity to obtain a wide range of valuable FIR 
fine-structure diagnostics never observed before in dwarf galax- 
ies. The Dwarf Galaxy Survey (DGS; P.I. Madden) is survey- 
ing the dust and gas of a wide range of low metallicity galaxies 
in the local univ erse dCormier et al.ll2010tlG alametz et al. 2010; 
lO'Halloran et al.ll2010tlMadden et al.l l2012b HRemv et al.ll2012l) 
using PACS (IPoglitsch et all 1201 Oh and SPIRE dGriffin & et al.l 
2010). The brightest low-metallicity galaxy of the DGS observed 
with the PACS spectrometer is Haro 11, an infrared-luminous 
starburst galaxy with metallicity Z = 1/3 Z . Haro 11 has the 
largest dataset obtained with PACS with observations of 7 FIR 
fine-structure lines: [C n] 157, [O in] 88, [O i] 63, [O i] 145, 
[N m] 57, [N n] 122, and [N n] 205 fim. It has also been ob- 
served with the IRS spectrograph onboard Spitzer. Unresolved in 
the IR, Haro 1 1 is an ideal candidate to conduct a study model- 
ing an exhaustive set of gas tracers to retreive information about 
the structure of the multi-phase low metallicity ISM. 

In this paper, we propose a systematic approach to model 
the multi-phase ISM of Haro 11, focusing on the gas proper- 
ties of the ionised and neutral ISM phases. We consider a set 
of 17 Herschel and Spitzer MIR and FIR ionic and atomic fine- 
structure cooling lines which trace a wide range of environmen- 
tal conditions, spanning excitation potenti als from 8 to 55 eV . 
With the spectral synthesis code Cloudy dFerland et al.l [l998). 
we create a self-consistent model of the dominant ISM compo- 
nents from which these lines originate in order to understand the 
prevailing mechanisms and physical conditions. In Section[2l we 
present spectroscopic data from the Spitzer and Herschel obser- 
vatories. Section [3] explains the modeling strategy and assumed 
initial parameters. Section [4] [5] and [6] report the model results 
and conditions of the ISM components. This study establishes 
an encompassing picture of the dominant ISM phases in Haro 1 1 
and describes the mass budget of the various phases in our model 
(Section|7]i. 

1.1. Studies of Haro 1 1 

Several studies of Haro 1 1 have highlighted its peculiarity as 
a galaxy and challenged our understanding of star formation. 
Haro 11 - also know as ESO 350-IG38 - is a blue compact 
dwarf (BCD) ga laxy at 84 Mpc, wi th M B = -20.3, and metal- 
licity Z~l/3 Z dGuseva et al.ll2012h . It is composed of 3 main 
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star- forming reg i ons, o r knots, described in IVaderetal J (1 19931) 
and[Kunthet al. (2003), w ith a mo rphology presumably result- 
ing from a merger event dOstlin et al.l 120011) . The starbursting 
nature of Haro 1 1 is visible in Figure Q] with a dominant bright 
stellar component and dust lanes in front of knot B. These knots 
host hundreds of young star clusters and ~60 super star clus- 
ters (SSC) with ages between 1 and 100 Myr pea king at 3 Myr, 
and c onstituting a total stellar mass of 10 10 M Q (lAdamo et al.l 
120101) . An older stellar population is also present, primarily in 
knots A and C, as seen in the red colors of the V-K bands, mod- 
eled in lMichevaetal.1 d2010l) by a metal-poor (Z=0.001) stellar 
population of age around 14 Gyr with standard Salpeter IMF, in- 
dicating that the gala xy is not undergoin g its first star formation 
event. Lya emission dKunth et alJI 19981 ) and Lyman continuum 
leakage were i nvestigated using UV and X-ray observat i ons in 
iBergvall et all (120061): paves et al] f2007[): iGrimes et alJ d2007l) 
and iLeitet et alj d2011l) . In iLeitet et al.l d2011l) . thev estimated 
the leakage to be 3%. In particular, the discrepant location of 
Lya and Ha emission, originating in the ionised gas surrounding 
the star-forming regions, may be explained by a diffuse ionised 
ISM in the halo of the galaxy with internal resonant scattering 
at the surface of c lumpy dense neutral clouds dKunth et alJ2003t 
lHaves et alj2007l) . These properties have resulted in its popular- 
ity as a local analogue o f the high-redshift Lyman-break galax- 
ies dOverzier et al . 2008), and therefore a must-study case to link 
nearby star formation to the distant universe and galaxy evolu- 
tion. 




Figure 1. HST/ACS 3-color image of Haro 1 1 with labels 
for the 3 star-forming knots (blue: F435, green: F550, red: 
FR656N/Ho'). Images were downloaded from the Hubble 
Legacy Archive. 



In the IR, Haro 1 1 is extremely bright compared to other 
nearby dwarf galaxies. IRAS measurements yielde d a FIR lu- 
minosity Lfir of 1 x 10 11 L P dSanders et al.ll2003ah from 40 to 
400 pm using the lLonsdale et alJ d 19851) prescription. Its spec- 
tral energy distribution (SED) shows a peculiar behavior. It 
peaks around 40 pm, indicating the prevalence of warm dust, 
and shows an excess of emission in the submm, as unveiled by 
AP£X/LABOCA and Herschel, whi ch is often seen in low metal- 
licity starbursting dwarf ga l axies dGalliano jtIDl2003ll2005l 
iGalametz et all I2009L l201lt iRemv et al.ll2012l) . The origin of 



the excess is not well det e rmine d and several hypotheses de- 
scribed in IGalametz et al.l (1201 ll) . including a change of dust 
emissivity, spinning dust grains, or a cold dust component, 
are i nvoked in current dus t models to explain the submm ex- 
cess. iGalametz et al.l d2009l) present a detailed dust modeling of 
Haro 1 1 in which they derive a total dust mass of 6 x 10 6 M if 
the submm excess is ignored, and of 2 x 10 7 M if the submm 
excess is interpreted as cold dust. Assuming the excess as cold 
dust (T~ 10 K) is however difficult to reconcile with the observed 
dust and gas mass budgets, and the exp ected lower dust-to-gas 
mass ratios (DGR; Gal ametz et alj|20lH) . Integrating the SED 
from 3 to 1100//m, th e total IR luminosity L TIR is 1.4 x 10 11 L 
(Galametz et al.l20"09l) . We will use this value of Ltir throughout 
this paper. 

Spectroscopic observations show bright lines in the FIR with 
/SO/LWS, while both the H 1 21 -cm line, and the CO (J=l-0) 
line, the mo st common tracer of molecular gas in galaxies, are 
undetected (Bergv all et al.ll2000T) . A detailed study of th e neu- 
tral gas of the ISM of Haro 1 1 by IBergvall et al] d2000l) inves- 
tigated the mass budget of the different ISM phases. They es- 
tablish masses for each phase: M(H 1) < 10 8 M , with an up- 
per limit on the atomic H 1; M(PDR) of 2 x 10 8 M using 
the LWS [C n] observations, the PDR mass being larger than 
the H 1 mass; and M(H2) ~ 10 8 M for the molecular gas, us- 
ing both a scaling of the PDR mass and the CO upper limit 
converted to H2 column density through the CO/H ? conversion 
factor (Xqo)- They used a Galactic value for Xqo (Stron get al.l 
fT988h : Z co , g ai = 2.3 x 10 20 cm 2 (Kkm s" 1 )" 1 . Therefore they 
describe the ISM structure as: very extended PDRs, in which the 
H 1 mass resides; a small fraction of the total ISM gas is left in 
atomic form; and the bulk of the gas, that is normally in ne utral 
state, is ionised by the starburst. Bergvall & Ostlin (2002|) de- 
rive M(H n) ~ 10 9 M from Ha observations for the ionised gas. 
These ma sses are lower than the total stellar mass of 10 10 M 
derived in lAdamo et all d2010l) . 

However, important questions remain, particularly on the 
state of the ionised, neutral atomic, and molecular gas. While 
evidence for young massive stellar populations and SSCs exist 
from the optical and IR point of view, evidence for any sub- 
stantial reservoir of neutral gas, in the form of either atomic or 
molecular, is non-existent to date. Where is the fuel for the vig- 
orous star formation observed in Haro 11? What is the state of 
the gas reservoirs in Haro 11, residing in the ionised, neutral, 
and molecular form? This study, using a wide variety of tracers 
of neutral and ionised gas, aims for a more complete description 
of the enigmatic nature of the ISM of Haro 1 1 . 



2. Observations 

2.1. Spitzer data 

Haro 1 1 was obser ved with the Infrared Spectrograph (IRS; 
iHouck et al.l 1 2 004) on board the Spitzer Space Telescope 
(IWerner et alE 004). on the 17th of July 2004 (PL: J. R. Houck). 
The 4 modules were used, i.e., the 2 high-resolution modules: 
Short-High SH (« 10-20 pm) and Long-High LH (« 19-37 /mi) 
with R ~ 600, and the 2 low-resolution modules: Short-Low 
SL (« 5 - 14 pm) and Long-Low LL (« 14 - 37 pm) with 
R ~ 60-120. Observations were done in staring mode (AORKey 
9007 104). Figure|2]shows the slits of the high-resolution module 
on the sky. 
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Table 1. General properties of Haro 1 1 . 



Quantity 



Value 



Reference 



RA (J2000) 
Dec (J2000) 
Distance 
12+log(0/H) 
Optical size 
i (<■) 

£fuv 

-^X-rays 



00h36'52.5" 
-33°33'19" 
84 Mpc 
8.2 

0.4'x0.5' 
1.4x10" L 
2 xlO 10 L G 
4.4xl0 7 L G 
1.9xl0 8 L o 
8xl0 8 L„ 



NED w 
NED 

NED (i) 

Guseva et al. (2012) 
NED 

Galametz et al. (2009) 
Grimes et al. (2007) 
Grimes et al. (2007) 
Haves et al. (2007) 
Ostlin etal. (1999) 



(a) NASA/IPAC Extragalactic Database, (b) This distance is derived 
with a Hubble constant Ho = 73 km/sec/Mpc. (c) Ltir is the integrated 
SED from 3 to 1 100 pm. All luminosities are scaled to the distance 
chosen in this paper of 84 Mpc. 




Figure 2. Ha image of Haro 1 1 from the Hubble Legacy 
Archive. The Spitzer/IRS high-resolution slits are overlaid in 
green and the Herschel/PACS array in red. The PACS beamsize 
at 150 pm is 1 1.5". The IRS low-resolution slits are not shown. 



The low-resolution spectra w ere retrieved from the CASSIS 
Atlas[3 (Lebouteiller et al. 2011), with a tapered spatial window 
that scales with wavelength to account for the broader PSF at 
longer wavelengths. The high-resolution spectra were obtained 
with SMART v8.2.2 (ILebouteiller et alJ l20Tfj iHigdon et all 
|2004|) from the full-slit extraction of SH and LH and assum- 
ing a point-like source calibration. The SL and SH spectra were 
stitched to LL and LH by applying the same factor 1.13, which 
is due to the fact that the source is slightly extended for the SL 
and SH apertures (3.7" and 4.7" height respectively), while it is 
point-like for the LL and LH apertures (« 11"). Another factor 
of 1.15 was applied to the low -resolution spectra to match the 
continuum of the high-resolution. The final rest-frame spectrum 
is displayed in Figure[3](fo/?), with a zoom on individual spectral 



1 The Cornell Atlas of Spitzer/IRS Sources (CASSIS) is a product of 
the Infrared Science Center at Cornell University, supported by NASA 
and JPL. 



lines (bottom). With a spectral resolution >500 km s _1 , the lines 
are not spectrally resolved. 

Line fitting and flux measurements were also done in 
SMART. A Gaussian and second order polynomial function 
were fitted to the line and continuum. For total uncertainties, we 
added 10% of the line flux due to calibration uncertainties. The 
resulting fluxes and uncertainties of the most prominent lines 
are listed in Table [2] The fluxes agree within the uncertainties to 
those of IWu et aijJ2006l) . The [Ne m] fine at 36.0 pm was not 
observed as it falls out of the IRS spectral range. 

2.2. Herschel data 

Haro 1 1 was obser ved with the PACS spectrometer 
dPoghtsch etal J 1201 Oh on the 27th of June 2010, OD 409 
of the Herschel mission dPilbratt et al.1 1201 Oh as part of the 
Guaranteed Time Key Program, the Dwarf Galaxy Survey (P.I.: 
S. C. Madden). The PACS array is composed of 5x5 spatial 
pixels each 9.4" square, covering a total field-of-view of 47". 
Haro 1 1 was mapped with 2x2 rasters separated by 4.5" in the 
7 following fine-structure spectral lines: [C n] 157 pm, [O in] 
88 pm, [O i] 145 pm (OBSID 1342199236), [O i] 63 pm, [N m] 
57 pm, [N n] 122 pm (OBSID 1342199237), and [N n] 205 pm 
(OBSID 1342199238), for a total of 7.2h. The [N n] 205 pm line 
was also re-observed on OD 942 for 1.6h with a single pointing 
(OBSID 1342234063). The observations were done in chop-nod 
mode with a chop throw of 6' off the source. The chopped 
position is free of emission. The beam size is 9.5" at 60 pm, 
and 11.5" at 150 pm (diffraction limited above ~100 pm). The 
spectral resolution is ~ 90 km s at 60 pm, 125 km s _1 at 
90 pm, 295 km s _1 at 120 pm, and 250 km s" 1 at 150 pm (PACS 
Observer's Manual 201 10). 

The data were reduced with the Herschel Interactive 
Processing Environment (HIPE) v7.0.1935 dOttll2010l) . We used 
standard scripts of the PACS spectrometer pipeline. The line fit- 
ting and map mak i ng we re done using the software PACSman 
(ILebouteiller et al.l 120121) . The signal from each spatial posi- 
tion of the PACS array is fitted with a second order polyno- 
mial plus Gaussian for the baseline and line. For total uncertain- 
ties, we add 30% of th e line flux due to calibration uncertainties 
(Poglitsch eTaT]l2010h . The rasters are combined by drizzling to 
produce final reduced maps with pixel size of 3". 

With Herschel, all of the lines are well detected, at least in 
the central pixel of the PACS array, except the [N n] 205 pm 
line. The [N n] 205 pm line is located at the edge of the first 
grating order of the spectrometer and affected by spectral leak- 
age, due to overlapping of the grating orders. It is contami- 
nated by signal from the second grating order around 100 pm, 
and is therefore difficult to observe and to calibrate. The line 
profiles are displayed in Figure [4] The line centers range be- 
tween -30 and +70 km s _1 within the maps of all observed 
lines, and indicate a rotation around the north-south axis of both 
the ionised and neutral gas components, slightly t ilted com- 
pared to the velocity structure of Ha analysed in Ostli n et all 
(1999). All lines display broad profiles, with a broadening up 
to 300 km s _1 , which suggests that both the neutral and ionised 
gas are coming from several components with different veloc- 
ities. However, with spectral resolution >100 km s , the lines 
are well fitted by a single Gaussian and we cannot separate 
the different velocity components in the spectra. This suggested 
complex velocity structure is corroborated by studies of opti- 
cal lines that show broad profiles as well dKunth et alJ fl998: 



available at http : //herschel . esac . esa . int/Docs/PACS/html/pacs_ 
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Figure 3. Top: Spitzer/IRS spectrum of Haro 1 1 in the rest wavelength frame (z = 0.0206), SL in green, SH in blue, and LH in 
red. Emission lines are labeled on top. The IRAC and MIPS photometry is also overlaid (fluxes in Table |3). Bottom: Zoom on the 
individual spectral lines. Fluxes and upper limits are reported in Table|2] 



lOstlin etal.ll999b , and where the neutral gas is shifte d compare d 
to the ionised ga s, expanding over few kpcs dKunfh et al.l l2003). 
FUSE spectra in Grime set al.l (120071) show absorption lines of 
the ISM with FWHM 300 km s' 1 and blueshifted by 80 km s 1 . 
The higher ionisatio n lines can be stronger in galactic outflows 
(IContursi et al.ll2012l) . 

The [C n] line shows a marginally extended emission at 
157 /mi, peaking at a signal-to-noise ratio of 50 in the center 
and dropping below detection at the edge of the map. The emis- 
sion of the other lines is rather compact. The shape of the PACS 
spectrometer beam is not known well enough to measure an ac- 
curate source size if the source were slightly extended. Hence we 



cannot recover spatial information. To derive total line intensi- 
ties, we integrated the fluxes over a circular aperture of diameter 
30" to encompass all the emission. The line intensities are listed 
in Table [2] Observations of the [O m], [C n], and [O i] 63 and 
145 /im lin es were also performed with the LWS instrument on- 
board ISO (iBergvall et al.ll2000t iBrauher et aDl2008b IVasta et al.l 
120101) . The PACS intensities of the [O i] 63 fxm and [O m] lines 
agree within 40% with the LWS observations, although they are 
found to be lower with PACS. This may be attributed to a beam 
effect, the LWS beam size being ~80". However, we find that 
the [C n] intensity is two times brighter with PACS compared 
to LWS. To understand this discrepancy, we have reduced the 
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Figure 4. Herschel/PACS line profiles from an individual spatial pixel (9.4") centered on the galaxy. 



PACS data with different flat-field corrections, as well as with 
background normalization (PACS Data Reduction Guide). All 
these methods give the same final [C n] flux as quoted in Table|2] 
within 10%, ruling out a mis-calibration from PACS. Therefore 
we attribute the difference between the PACS and LWS fluxes to 
a calibration problem from the LWS instrument. 



2.3. Spectral line description 

The IRS and PACS spectra probe the ionised and neutral media. 
The IRS spectrum shows a wealth of MIR fine-structure lines, 
the brightest being (in order) [Ne in] 15.6 yum, [Si n] 34.8 //m, 
[S in] 18.7 and 33.5 /im, [S iv] 10.5 /mi, and [Ne n] 12.8 /mi. 
PAH features at 6.2, 7.7, 8.6, 1 1.2, and 12.8 /jm are also present. 
The PACS spectrum also shows bright FIR lines, especially 
[O m], [C n], [O i] 63 /mi, and [N m]. The [O m] line is the 
brightest of all, with [O m]/[C n] = 2.4. While [C n] is the 
brightest FIR fine- structure li ne in normal and dusty starburst 
galaxies (Bra uher et al.l 120081) . we find th at this is often not 
the case in low-me t allicity dwarf galaxies (|Hunter et al] |200 1 ; 
ICormier et aIll2010ll2012aHMadden et alJl2012bh ~ 

The young star-forming nature of this galaxy also agrees 
with intense [Ne m] (41.0 eV), [O m] (35.1 eV), [S iv] (34.8 eV), 
and [N m] (29.6 eV) lines observed, arising from relatively 
highly ionised regions. [O iv] (54.9 eV) at 25.9 /urn, which re- 
quires a relatively hard radiation field (Sect. 16.31 . is also de- 
tected. The [Ne m]/[Ne n] ratio is 3.5, [S iv]/[S m] is 0.9, and 
[N m]/[N n] is 8. These r atios are indicative of a hard radia- 
tion fi eld typical of BCDs dThornlev et all 120001: iMadden et al] 
2006). While [C n] can come from both the diffuse ionised phase 
and the neutral phase, [O i] traces only the neutral phase. The 
[O i] 63 /mi/[C n] ratio is 0.9, s imilar to that found in other dwarf 
galaxies and dustier galaxies dHunter et al.ll2.Q0ll : iBrauher et alJ 



2008). Table|2]lists the ionisation energies of the observed lines, 
as well as the excitation temperatures and critical densities (n c „>) 
of the transitions with collisions with H and/or e~ . 

Comparing line luminosities to the TIR luminosity, we find 
L [C iij/Ltir = 0. 1 %, L [0 mi/LuR = 0.3%, and L IO rj/Lrm = 0. 1 %. 
Altogether, the FIR lines represent 0.6% of the TIR luminosity; 
and when adding the MIR lines, they represent 1.2% of Ljir. 
Together these lines are responsible for most of the gas cooling 
of the IS M. Such line-to - TIR r a tios are typical of st ar-forming 
galaxies dMalhotra et all 120011: IBrauher et all [200810. Similar 
Lrc iii/Ltir ratio s are found in several h igh-redshift g alaxies 
(Mai olino et all 120091: IStacev et alJbOlOh . ICroxall et al.l d2012l) 
find L[c ii]/L T ir ~ 0.5% and Lp i]/L T ir ~ 0.1% in spatially re- 
solved regions of the nearby spiral galaxy NGC 4559. There 
is observational evidence that L[c iij/Ltir tends to decrease in 
galaxies with warmer FIR-color, f v (60 //m)/f v (100 /mi). Among 
other effects, this may be attributed to a less efficient photo- 
electric heating du e to charged grains (Malh otra et al.l l2001 : 
ICroxall et al.l 1201 % , or to a UV-photon screening from large 
amounts of dust as in ultra-luminous IR g alaxies (ULIRGs) 
dAbel et all l2009t iGracia-Carpio et al.l l201ll) . With observed 
FIR-color ~1, Haro 1 1 falls within the large spread of L[c ii]/Ltir 
values and is not particularly [C n] defici ent nor [C n] excessiv e, 
as often seen in low metallicity galaxies (Madden et al. 2012b). 

2.4. Photometry data 

For the discussion on the SED of Haro 1 1 in Section lTTl we use 
broad-band data available from the X-ray to the radio regime, 
to assess how the models perform in reproducing the contin- 



3 Note that they u se L FIR , determined from IRAS fluxes with 
the iHelou et al.l fl988h formula, and that Ltir estimated as 2 x Lfir 
dHunter et alj|2001l) can be different from the true L T i R when including 
submm observations. 
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Line 


Wavelength 


Flux (fl) 


Uncertainty 


Excit. potential^ 


Excit. temperature^ 1 


Critical density*''' 




Gum) 


(10~ 16 W m" 2 ) 


(10~ 16 W nr 2 ) 


(eV) 


(K) 


(cm" 3 ) 


Spitzer/IRS 


[Siv] 


10.51 


5.37 


0.66 


34.79 


1 369 


5 xlO 4 [e] 


[Ne n] 


12.81 


3.53 


0.43 


21.56 


1 123 


7 xlO 5 [e] 


[Ne m] 


15.56 


12.3 


1.5 


40.96 


925 


3 xlO 5 [e] 


[S m] 


18.71 


5.79 


0.72 


23.34 


769 


2 xlO 4 [e] 


[S m] 


33.48 


5.85 


0.69 


" 


430 


7 xlO 3 [e] 


[Si ii] 


34.82 


6.21 


0.75 


8.15 


413 


3 xlO 5 [H], 1 xlO 3 [e] 


Humphreys a 


12.37 


0.208 


0.045 


13.60 


1 163 




H2 (0,0) S(0) 


28.22 


<0.62 




4.48 


510 


7 xlO 2 [H] 


H2 (0,0) S(l) 


17.03 


0.168 


0.055 




1015 


2 xlO 4 [H] 


H2 (0,0) S(2) 


12.28 


0.101 


0.013 


» 


1682 


2 xlO 5 [H] 


H2 (0,0) S(3) 


9.66 


0.121 


0.018 




2504 


9 xlO 5 [H] 


[Fe m] 


22.93 


0.54 


0.20 


16.19 


627 


1 xlO 5 [e] 


[O iv] 


25.89 


0.49 


0.18 


54.94 


555 


1 xlO 4 [e] 


[Fe n] 


25.99 


<0.25 




7.90 


554 


2 xlO 6 [H], 1 xlO 4 [e] 


[Ar m] 


8.99 


0.99 


0.33 


27.63 


2 060 


3 xlO 5 [e] 


[Aril] 


6.99 


0.60 


0.15 


15.76 


1 600 


4 xlO 5 [e] 


[Ne v] 


14.32 


< 0.080 




97.12 


592 


3 xlO 4 [e] 


HerschelfPACS 


[C ii] 


157.74 


7.4 


2.2 


11.26 


91 


3 xlO 3 [H], 50 [e] 


[O m] 


88.36 


18.0 


5.4 


35.12 


163 


5 xlO 2 [e] 


[O i] 


63.18 


6.7 


2.0 




228 


5 xlO 5 [H] 


[O i] 


145.52 


0.56 


0.17 




99 


1 xlO 5 [H] 


[N m] 


57.32 


2.50 


0.78 


29.60 


251 


3 xlO 3 [e] 


[Nil] 


121.9 


0.31 


0.10 


14.53 


118 


3 xlO 2 [e] 


[Nil] 


205.18 


< 0.29 






70 


45 [e] 



The /50/LWS fluxes can be found in Beravall et al. (2000), Brauher et al. (2008), and Vasta et al. (2010). 

(a) Upper limits are given as the 3-<r uncertainty from the line fit. 

(b) Excitation potential. Energy to create an ion. For Hi, it is the energy to photodissociate the molecule. 

(c) Excitation temperature. Temperature T=AE/k required to populate the transition level. 

(d) Critical den sity n erlf noted [H] for c ollisions with hydro gen a toms (T=100 K and T =300 K for the H2 lines), [e] with electrons (T=10 000 K). 
Values are from Malh otra et alj ( l20011 ). lGiveon et alj ( 120021) . and I Kaufman et all (12006). 



uum emission. The photometry is taken from N ED. In addition , 
we use the Spitzer/MlPS photometry from Ben do etaD (l2012[). 
the H erschel/PACS and SPIRE photometry from iRemv et al] 
d2012l) , a nd the SpitzerflRAC an d LABOCA 870 pm photom- 
etry from Gala metz et al.l (120091) . All fluxes are summarized in 
Tabled 



3. Modeling strategy 

3.1. Identification of the model components 

We aim to model the IR emission from all 17 detected lines (ex- 
cept the MIR H2 lines which will be modeled in a coming paper) 
listed in Table |2]to understand the conditions of the various me- 
dia and the dominant processes at work in this low-metallicity 
starburst. The ISM is very heterogeneous, and over the scale of a 
galaxy, various phases with a wide range of conditions are mixed 
together. 

In our spectral dataset, we notice that Haro 1 1 is typified by 
very intense [Ne m], [Ne n], [S 111], [S iv], [O in], and [N in] lines. 
These high-excitation lines require a hard radiation field which 
is naturally encountered in a compact H n region. Stars in com- 
pact H 11 regions are young, hence their UV radiation is harder, 
and SSCs are likely to contain massive stars which contribute to 
the production of hard photons. Moreover, the high critical den- 
sities of these lines (except the [O m] line) favor their origin at 



the same location. Therefore, the intensity of these lines demon- 
strates the importance of the dense H n region in the global 
spectrum of Haro 1 1 . Moreover, the Herschel observations show 
bright [O 1] and [C 11] lines, which are usual tracers of the neu- 
tral gas. In particular, the intensity of the [0 1] lines indicates the 
presence of dense PDRs. These two dense (ionised and neutral) 
phases will be investigated first since they account for the emis- 
sion of most of the spectral lines. Derived physical conditions 
are presented in Section[4] 

With low critical density (~300 cm 4 ) and low ionisation 
potential (14.5 eV), the [N 11] 122 pm line is a tracer of low- 
excitation diffuse gas rather than the dense H 11 region. The in- 
tensity of [N 11] demonstrates the need for an additional diffuse 
low-ionisation component that is investigated in Section [5] Our 
subsequent modeling also shows that this diffuse phase is the 
main contributor to the [Ne 11] and [C 11] lines. 

Therefore, the three main phases that we identify as neces- 
sary model components to account for the observed emission 
lines are: 1) a dense ionised phase (Sect. 14. IT ). 2) PDR (Sect. l4.2l l. 
3) a diffuse low-ionisation medium (Sect. [5]). 

In addition, we discuss the possible presence of other model 
components that may influence the line predictions. The detec- 
tion of the high-excitation [O iv] line is the driver of a discus- 
sion on the role of X-rays and XDRs in Section [63] The [O m] 
line, which is the brighest of all MIR-FIR lines, has high exci- 
tation potential (35.1 eV) and low critical density (~500 cm 4 ). 
Although it is accounted for by the dense H 11 model, it may orig- 
inate as well from diffuse high-excitation medium that we do not 
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Table 3. Photometry data. 



Instrument 


Wavelength 


Flux 


Uncertainty 


Referer 






(Jy) 


(Jy) 




Chandra 


2-10keV 


9.17e-09 


1.03e-09 


(1) 


Chandra 


0.5-2 keV 


3.21e-08 


0.43e-08 


(1) 


ROSAT 


0.1-2 keV 


4.41e-08 


0.55e-08 


(2) 


FUSE 


1150 A 


2.30e-03 


0.33e-03 


(3) 


GALEX 


1530 A 


3.63e-03 


7e-05 


(4) 


GALEX 


2315 A 


4.37e-03 


4e-05 


(4) 


Cousins 


0.439 fim 


8.03e-03 


0.69e-03 


(5) 


Cousins 


0.639 yum 


0.012 


0.001 


(5) 


2MASS 


1.25 yum 


0.013 


0.001 


(6) 


2MASS 


1.64 yum 


0.013 


0.001 


(6) 


2MASS 


2.17 fim 


0.014 


0.001 


(6) 


IRAC 


3.6 /jm 


0.023 


0.002 


(7) 


IRAC 


4.5 yum 


0.029 


0.003 


(7) 


IRAC 


5.8 fim 


0.073 


0.007 


(7) 


IRAC 


8 yum 


0.177 


0.018 


(7) 


IRAS 


12 fim 


0.52 


0.02 


(8) 


MIPS 


24 ftm 


2.51 


0.03 


(9) 


IRAS 


25 fim 


2.49 


0.03 


(8) 


IRAS 


60 fim 


6.88 


0.04 


(8) 


MIPS 


70 fim 


4.91 


0.49 


(9) 


PACS 


70 fim 


6.08 


0.18 


(10) 


IRAS 


100 yum 


5.04 


0.03 


(8) 


PACS 


100 fim 


4.97 


0.15 


(10) 


MIPS 


160 yum 


2.01 


0.24 


(9) 


PACS 


160 yum 


2.43 


0.12 


(10) 


SPIRE 


250 yum 


0.633 


0.042 


(10) 


SPIRE 


350 yum 


0.231 


0.013 


(10) 


SPIRE 


500 yum 


0.092 


0.008 


(10) 


LAB OCA 


870 yum 


0.040 


0.006 


(7) 



(1) Haves et al. (2007); (2) Taier et al. (2005); (3) Grimes et al. (2009) 
(4) [jalesias-Paramo et al. ( 200 6}): (5) | Lauberts & Valentiinl 
(6) Uarrett et alj d200d)~77) iGalametz et al 
d2003bl) : (9) lBendo et all d2012h : QOI lRemv et alJ d2012h . 



(1989); 

J2009I); (8) ISanders et al] 



model but briefly discuss in Section 1672] Finally, we investigate 
the influence of magnetic fields on the cloud density profile and 
their impact on the PDR predictions in Section l6~T1 



to the inner edge of the cloud is the inner radius (r; n ). Depending 
on where the calculation is stopped, the ionised, atomic, and 
molecular phases can be treated as the model is computed within 
one Cloudy model. Here we stop the models until the molecular 
fraction H2/H reaches 50% such that both the ionised and PDR 
phases are computed within the same model. We opt to impose 
pressure equilibrium throughout the models. This is particularly 
important for the density profile between different phases. We 
will come back to this assumption in Sect. 14.21 



3.2.1 . Parameter grid of the dense H 11 region 

We vary the following 3 physical parameters in order to calculate 
our model grid: the age of the burst (Abmst), the initial hydrogen 
density at the inner edge of the cloud (ne), and ri n , where the 
calculation starts (see Fig|5]for an illustration). The age is varied 
from 2.4 to 5.0 Myr (Se ct. 143721 r in from 10 200 to 10 22() cm 
(or 3 pc to 3 kpc, Sect. 147131 . and n H from 10 1 to 10 4 cm 3 
(Sect. l47TTn >. The step size for Ab ms t is 0.2 Myr (linear step), and 
for ri n and nn, we use 0.2 step sizes on a logarithmic scale. 






min 


max 


step 


unit 




2.4 


5.0 


0.2 


Myr 




1 


4 


0.2 


log cm 3 


r in 


20 


22 


0.2 


log cm 



Figure 5. Geometry of the Cloudy model with the three free 
parameter that we vary in the grid: the age of the burst (Aburst), 
inner radius (rj n ), and initial hydrogen density (ne). Parameter 
boundaries and step sizes are also indicated. 



3.2. Method 

We use th e ID spectral synth esis code Cloudy vlO.00, last de- 
scribed bv lFerland et alj d!998l) . which includes p hotoionisation 
and photodissociation treatment (lAbel et al.ll2.005h . to model the 
multiple phases from which the emission lines originate. Since 
Haro 1 1 is not resolved in most of the Spitzer and Herschel 
bands, we have no spatial constraints. Thus, the model is re- 
duced to a central source of energy surrounded by a spherical 
cloud. In the case of Haro 11, this central source is a starburst. 
The dense ionised gas that is closest to the central source with 
respect to the other phases that we model separately will yield 
the best constraints on the properties of the stellar cluster that is 
powering the model. Therefore we first derive the best fit model 
concentrating on the dense H 11 region diagnostics, which are the 
[Ne n] 12.8 yum, [Ne m] 15.6 /urn, [S rv] 10.5 fim, [S m] 18.7 yt/m 
and 33.5 yum, and the [N in] 57 yum lines. Moreover, the interpre- 
tation of these ionic lines is less ambiguous than, for example, 
the [C n] line. Following this step, we will assess what this model 
predicts for the PDR phase. 

The starburst ionises the inner edge of the cloud, where the 
H n region begins, and the radiative transfer is computed step by 
step progressively into the cloud. The distance from the source 



3.2.2. x 2 evaluation 

In order to find a solution for this ionised phase, we aim first to 
reproduce as best as possible the intensities of the ionic lines. 
For each individual model, the goodness of the fit is evaluated 
by computing the reduced x 1 - The x 1 is measured from a set 
of selected lines. The lines we consider to constrain the fit are 
the following: [Ne 11] 12.8 fim, [Ne m] 15.6 fim, [S iv] 10.5 fim, 
[S in] 18.7 yum, and 33.5 yum, because they are the brightest and 
relatively well understood, as well as the [N m] 57 fim line from 
PACS which is also bright and expected to arise from the same 
medium. With 6 observations and 3 free parameters to fit, the de- 
gree of freedom is N = n D b s - n param = 3. Best fit models are de- 
termined by minimising the reduced^ 2 . There are several lines 
that originate in part from the dense ionised gas that we do not 
consider when calculating the ^ 2 , although we do discuss how 
the model predictions of these lines compare to the observations. 
We do not consider the high-excitation [O iv] line because it is 
faint and can be excited by several mechanisms (see Sect. 16.3b . 
We do not include the [O in] and [N 11] lines here as constraints 
for the dense H 11 region because our subsequent modeling shows 
that a fraction of these lines arises from a low density ionised 
medium which we consider in a later step (see Sect. [5j. 
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3.3. Setting the input fixed parameters 

3.3.1. Incident radiation field 

We use Starburst99 v6.0.2 (iLeitherer et al.ll2010l) to reproduce 
the incident, wavelength-dependent radiati on field density of 
Har o 1 1 . Following the U V-optical studies bv lHaves et al.l 12007) 
and lGrimes et al.l (|2007), we assume a Salpeter IMF from 0.1 to 
100 M and stellar trac ks from Padova AG B with metallicity 
Z = 0.004. According to lAdamo et al.l d2010l) . the stellar popu- 
lation of Haro 1 1 is very young. The present starburst started less 
than 40 Myr ago, and shows a formation peak at about 3.5 Myr. 
Knot B is dominated by a young starburst of 3.5 Myr, and knot C 
contains an older starburst of 9.5 Myr (Fig. [I). Therefore, we se- 
lect the stellar spectrum of an instantaneous burst of a few Myr, 
letting the age of the burst be a free parameter, varying from 2.4 
to 5.0 Myr. We also add to this the spectrum of an instantaneous 
burst of 9.5 Myr to account for the slightly older stellar popula- 
tion. We keep this age fixed since the ionising part of the stellar 
spectrum does not vary significantly for ages larger than 5 Myr, 
and will therefore not significantly affect the model predictions 
of the ionised gas. This starburst spectrum dictates the shape of 
the input energy source in Cloudy. The incident radiation field on 
the cloud is shown in purple in FigureQ~2] For the input luminos- 
ity of the central source, we take the total observed infrared lu- 
minosity of Haro 11: Ltir = 1.4 x 10 11 L . This approximation 
assumes an energy balance between the UV and IR domains, i.e. 
that all of the UV photons are reprocessed in the IR by the gas 
and the dust. The emergent SED (bottom panel of Fig.[T2l shows 
that this approximation is reasonable since only 20% of the total 
emergent luminosity is emitted in the UV-optical, and 80% in the 
IR. We discuss in more details the energy balance in Sect l7.ll 

Haro 11 is detected in X-rays (e.g. iGrimes et alJ 

120071) . Previous studies on X-ray observations indicate 
the presence of soft (Lo.5-2kev = 1.9 x 10 8 L ) and hard 
(L2-iokev = 2.6 x 10 8 L Q ) X-rays. Therefore we add to the input 
SED an X-ray component that reproduces the shape of the ob- 
served X-ray spectrum and of total luminosity 5 x 10 8 L . The 
effects of the X-rays in the models are discussed in Section ROl 

3.3.2. Elemental abundances 

The elemental abundances directly impact on the predicted 
line intensities. To match the observed abundances in Haro 1 1 
as closely as possible, we use the abundances in the litera- 
ture. Abundances of oxygen, nitrogen, neon, sul phur, and ar- 
gon w ere obtained in knots B and C of Haro 1 1 bv lGuseva et al.l 
d2012l) using VL77X-shooter observations in the optical. In par- 
ticular, they find an oxygen abundance 12+log(0/H) of 8.1 in 
knot B and 8.3 in knot C, giving a O/H twi ce as large as the 
value 12+log(Q/H)-7.9 previously found by Bergvall & Ostlin 
(2002), altering the determination o f the metallici t y from ~l/6 
to 1/3 Z . We use the abundances in lGuseva et al.l (12012 ). aver- 
aged over the two knots, which have similar luminosities, for 
our calculations. Uncer tainties on these ab undances are 10%. 
For neon and sulphur. | Guseva et alJ (1201 2l) find similar val- 
ues to those found in IWu et al.l (120081) . derived from Spitzer 
data. The relatively hi gh neon abundance is typical of BCDs 
dlzotov & Thuanl["l999l) . For nitrogen, such high abundance is 
observed in a few BCDs (e. g . iThuan et al.ll996tlIzotov & Thuanl 
119981: iPustilnik et al.ll2004l). and can be explained b y the pres- 
ence of Wolf-Rayet stars dBergvall & Ostlinl 12002). which in- 
ject matter enriched in N and C into the ISM. For the other 
elements, we use gas-phase relative abundances typical of 



H ii regions, based on an average of the abundances in the 
Orion Nebula determined bvlBaldwin et all (1 19911) : iRubin etafl 
(ll991l) : IOsterbrock et all dl 9921) : lRubin et al.l(fl993l) . and scaling 
them to the metallicity of Haro 1 1 . The resulting carbon abun- 
da nce is log(C/Q) = -0.1 5, which is relatively high compared 
to llzotov & Thuanl d!999l) . where log(C/0) ~ -0.5. However, 
Bergvall et al. (2006) find that the ratio of C and O column den- 
sities is 0.3 (log) from the analysis of absorption lines in IUE 

and FUSE data. Such high value may indicate an enrichment 
in C and agrees with our adopted carbon abundance (assuming 
Orion-like C/O) within the uncertainties. For silicon, its abun- 
dance Si/O is low in the Orion abundance set compared to that 
usually found in BCDs, th erefore we use the aver age abundance 
of log(Si/0) = -1.5 from llzotov & Thuanl (fl999h . We consider 
the uncertainty on these theoretical abundances to be of a factor 
of 2. Table |4] summarizes the elemental abundances (and uncer- 
tainties) we use in Cloudy. The lack of measured abundances for 
the other elements limits their use to constrain the models (e.g. 
Fe, Si). 



Table 4. Elemental abundances set in Cloudy. 



Gas-phase composition 


Abundance (uncertainty) in log 


O/H 


-3.79 (0.14) 


NTH 


-4.64(0.12) 


Ne/H 


-4.44 (0.07) 


S/H 


-5.57 (0.10) 


Ar/H 


-6.20 (0.13) 


C/H 


-3.92 (0.3) 


Si/H 


-5.80 (0.3) 


Fe/H 


-5.92 (0.3) 


Grain composition 


Mass fraction 


Carbonaceous 


36% 


Silicate 


60% 


PAH 


4% 


Dust-to-Gas mass ratio 


2.64 x IQ- 3 



Abundances of oxygen, nitrogen, neon, sulphur, and argon are from 
iGuseva et ID ( 120121) . silicon from llzotov & Thuart d!999h . and the 
others are average abundances of the Ori o n Nebula determined by 
Baldwin et all d 1991h ; [Rubin et all fl99ll) : IOsterbrock et all dl992h : 
Rubin et alJ |l993:)> scaled to the metallicity of Haro 1 1 . 
The main elemental composition 0:C:Si:Fe in the grains is 4:9:1:1. 



3.3.3. Dust properties 

It is important to accurately set the dust properties for two 
reasons: (1) dust plays an important role in heating the gas 
through the photoelectric effect, which depends on the dust 
abundance and grain size distribution, and (2) the radiative 
transfer is strongly influenced by the properties of the grains. 
We first tried using the standard MRN grain size distribution 
(Mathis, Rump l. & Nordsieckll 1977b . which failed to reproduce 
the observed dust emission in the MIR, falling below the pho- 
tometry data by more than an order of magnitude. The MRN 
distribution underpredicts the number of small grains (minimum 
grain radius of 5 nm) and is usually not appropriate for mod- 
el ing the MIR emissio n in galaxies. Standard dust models such 
as [Draine &Li (2001) and lZubkoet al. I d200l do include more 
small grains than th e standard MRN model. Dust modeling of 
dwarf galaxies (e.g. iGalliano et al.ll2003l 120051 : iGalametz et al.l 
1201 lb also demonstrates the need for an increased abundance of 
smaller grains. Therefore, in Cloudy, we use the MRN distri- 
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bution extended down to 1 . 1 nm to account for the very small 
grains. Other grain size distributions have been used to model 
the IR emission of galaxie s, such as grain shattering distri- 
butions ( Donit aet al.l 120051) . that we discuss in S ection 17.1.21 
We also use PAHs as described in lAbel et all d2008l) . Both 
grains and PAH abundances are scaled to the metallicity of 
Haro 11 (1/3 Z ), fixing the dust-to-gas mass ratio (D/G) to 
2.64 x 10~ 3 (-1/400, the Galactic D/G is -1/150). This yields a 
ratio of visual extinction to hydrogen column density, Ay/N(H), 
of 10~ 22 mag cirT 2 . 



3.4. Limitations of the approach 

Our first limitation lies on the geometry of our modeling that 
we aim to constrain as much as possible with diagnostics from a 
large set of spectral lines. Some of the derived parameters of the 
compact H n region are strongly dependent on the (simple) ge- 
ometry we impose on the model, which we will see is a recurring 
theme throughout the paper. The parameters that are best con- 
strained, irrespective of the geometry, are the electron density, 
temperature and mass of ionised gas. As found in Section l4.1.3l 
the geometrically thin ionised gas shell, in combination with the 
large inner radius (~1 kpc) that describes the model solution of 
the dense H n region, seems physically implausible. Despite the 
geometry of the model, we refer to it with the term compact 
since this model component represents what is observationally 
described as a compact H n region. This is a direct result of as- 
suming a single region with a single central source. Although we 
are aware that the knots B and C of Haro 1 1 are spatially sepa- 
rated in the optical, we do not have the spatial resolution in the 
Spitzer and Herschel observations to treat them separately. If, for 
example, we were to construct a model as the sum of several re- 
gions - as seems reasonable when perusing the optical images - 
the inner radius of each region would have to decrease in order to 
obtain a similar effective U parameter, while the shell thickness 
would increase in order to obtain the same gas mass, thus evolv- 
ing towards a model with a more geometrically thick layout. It is 
impossible to obtain reliable and reproducible results using such 
a complex model, without better observational constraints on the 
actual source distribution, i.e. spatially resolved spectroscopy. In 
this paper we focus on the derived parameters that are least sensi- 
tive on the actual assumed geometry. In Section ITT1 we discuss 
the spectral energy distribution resulting from the combination 
of models of the various phases in terms of energy balance, but 
we do not pretend to accurately fit the total SED because of these 
geometry issues. 

Our second limitation lies in the evalutation of the goodness 
of the fits. It is mainly weighted by the elemental abundance 
uncertainties, which can be up to a factor of 2, and by the obser- 
vational uncertainties (up to 40%). Therefore, when comparing 
the model predictions to observed intensities, we aim for a match 
within a factor of 2. 



4. A model of the H n region and PDR: % 

4.1. The compact H n region 

The results of the computed grid are displayed in Figure [6] It 
shows the intensity predictions of the compact model labeled ^ 
for each ionic line as a function of density, ne, and inner radius, 
ri n , for a starburst age of 3.7 Myr (age of the best fit model, see 
Sect. l4.1.2l >. In this section, we discuss the influence of each grid 
parameter on the model predicions and possible solution values. 
These parameters are linked via the quantity U, which is the ion- 



isation parameter, defined as the dimensionless ratio of the inci- 
dent ionising photon surface density to the hydrogen density n^ 
U — 4^2^ ~ > where Q(H) is the rate of hydrogen-ionising pho- 
tons striking the illuminated face and c is the velocity of light. U 
describes how many ionising photons arrive per atom. 
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Figure 6. Predicted intensities for model ^ (l mo dei) normalised 
by their observed value (l b s ) for all ionic lines considered, as a 
function of ne and ri n for a given starburst age (of 3.7 Myr). 
The vertical dashed line shows the best fit model density 
ne =10 2 ' 8 cm 4 , and the vertical dotted line the critical density 
of the ionic specie, if within the range displayed. The best fit 
inner radius is r; n =10 2L4 cm (green symbols). 
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4.1 .1 . The density of the gas 

In general terms the predicted line intensities are dependent on 
the gas density through both the critical density of each atom, 
and the ionisation parameter, U. Line intensities increase with 
increasing gas density until the critical density of a fine-structure 
line is reached. Above such critical density, the de-excitation 
process is dominated by collisions rather than radiative decay. 
Therefore the maximum gas density is well determined by the 
critical density of the brightest line tracers. The gas density also 
plays an important role in determining the population of the dif- 
ferent ionisation stages. The ionisation parameter, U, is inversely 
proportional to the density (bottom right panel of Fig. [5}. 

We find that the minimum x 2 occurs for a solution with 
a density around 10 28 cirT 3 (dashed vertical line in Fig. 
This derived gas density compares well with density diagnos- 
tics proposed in the literature, for example , based on the rati o 
of the [S in] fines at 18.7 and 33.5 /im (iHouck et al l 11984 . 
This ratio is sensitive to the density since these lines origi- 
nate from the same ion but with significantly different criti- 
cal densities (Table |2j. From IXbel et al.l (|2005) and requiring 
[s iii] 33 ~ m ~ 1-0' we obtain a density of the ionised gas of 
nn = 10 27 cm 3 . The valu e of the ele c tron d ensity predicted 
by the diagnostic plot in IHouck et al.l d 19841) is also around 
10 28 cm" 3 . Following the method from iKingsburgh & Barlowl 
d 199211 19941) who use the programs EQUIB and RAT 10 (devel- 
oped at UCL by I. D. Howarth and modified by S. Adams), we 
also find nn ~ 10 2 8 cm -3 . 



4.1 .2. The age of the burst 

The spectrum of the starburst (intensity and shape) changes dras- 
tically with time in the age range considered here. The age of the 
burst sets the hardness of the incident radiation field. Figure [7] 
shows line ratios of [Ne m]/[Ne n], [O iv]/[Ne m], [O ra]/[N m], 
[O m]/[S iv], and [S iv]/[S m] 18.7 fim as a function of burst age. 
The observed ratios agree with an age of burst of 3.5-5 Myr. 
These values are also in agreement with the age of 3.5 Myr 
for the brightest knot B and 9 Myr for k not C, and a mixtur e 
of very young and older bursts, found in [Adamo et alJ ([2010). 
After 5 Myr, the ionising part of the burst spectrum decreases, 
the high ionisation lines start to disappear, and in particular the 
[Ne m]/[Ne n] ratio drops off (bottom panel of Fig. 13). This is 
why we do not compute models with Ab urs t > 5 Myr in our grid. 
The best fit on the neon lines favours a relatively young burst, 
of about 3 Myr. However, the emission of the highest excita- 
tion potential [O iv] line (54.9 eV) is rather favoured by models 
with ages between 3.7 and 4.7 Myr. Indeed this time frame cor- 
responds to the time needed for the onset of Wolf-Rayet stars. A 
burst age of 4.7 Myr is also a good match to the [Ne m]/[Ne n] 
ratio, although it does not agree with the intensity of the [S iv] 
line, which is one of the lines that has the highest weight in the 
X 1 - The best single solution model in that case has an age of burst 
of 3.7 Myr. 



4.1 .3. The inner radius 

All lines are sensitive to the inner radius r; n . Its choice has no di- 
rect influence on the hardness of the radiation field, but rather on 
its intensity. With increasing radius, the total input luminosity is 
distributed over a larger surface, decreasing the number of ion- 
ising photons per surface area, and therefore the ionisation rate. 
U decreases as the inverse square of the radius. Intensities of 
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Figure 7. Model predicted intensity ratios as a function of burst 
age Aburst, for a hydrogen density of 10 2 8 cm" 3 , and inner radius 
of 10 2L4 cm. The true burst age of these models is an average of 
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the vertical dotted lines the age of 3 Myr and 4.7 Myr that would 
better match the [Ne m]/[Ne n] ratio. 
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species with a high ionisation potential, such as [O iv], decrease 
with increasing r; n , and the inverse is observed for low ionisation 
species such as [N n]. For r; n < 10 210 cm, atoms are preferably 
present in high ionisation states. Models with such small ri n of- 
ten predict too much [O iv] and [S rv], and not enough [Ne n] 
and [S m]. [Ne m] and [S iv] intensities are constant at low rj n 
and then drop with increasing r; n . [S m], [N m] and [O m] inten- 
sities increase with low r; n (< 10 2L5 cm) because of an ionisation 
balance between S 2+ , N 2+ , 2+ and S 3+ , N 3+ , 3+ . Then their in- 
tensities drop with increasing rj n . The opposite effect is observed 
for the lower ionisation state species such as [Ne n] and [N n] 
which increase with r; n . 

The observations are best matched for an inner radius of 
10 214 cm (800 pc). At that distance, the model cloud is a very 
thin shell of large radius, at the center of which resides the cen- 
tral starburst. The resulting geometry is plane-parallel. This is 
not a realistic geometry, as discussed in Sect. 13.41 but we do not 
aim to model the spatial structures in Haro 1 1 . 



4.1 .4. Parameters of the best fit model 



'HII 



We have minimised the x 2 to find the parameters of the best fit 
model for the compact H n region described above. Figure [8] 
shows contour plots of the reduced x 2 computed for our grid 
of models. Contour levels at 1-cr, 3-cr, and 5-cr from the mini- 
mum x 1 are displayed in orange. The intensities of the spectral 
lines included in the fit are best reproduced with the parameters: 
n H = 10 2 8 cm" 3 , Abum = 3.7 Myr, and r in = 10 2L4 cm. The ioni- 
sation parameter is U = 10~ 2 ' 5 . The reduced^ 2 of this model so- 
lution, which has the minimum^ 2 , is 1.48. This value is close to 

1, indicating that this model describes the data well, with a prob- 
ability of ~70%. For clarity throughout the following Sections, 
we will refer to this model solution of the compact H n region as 

&H1I- 

With the strongest MIR line, [Ne m], is reproduced 

within 15%, as well as the [O iv] line. The two [S m] lines 
are achieved within 5% and their ratio within 1%. The modeled 
[S iv] and [Ar m] lines deviate by 20%. The [N m] line is repro- 
duced within 40%. The [Ne n] line, however, is under-predicted 
by a factor of 5 by 

We then compare the line intensities predicted by ^hii to the 
PACS lines, [O m] and [N n] 122 /mi, which were not used to 
constrain the model fit since they can originate from less dense 
ionised gas. We find that [O m] is under-predicted by a factor of 

2, which is within our expectations, while the [N n] 122 /mi is 
under-predicted by a factor of 4. The dense H n region is not the 
primary source of [N n]. The prediction for the [N n] 205 fim 
line falls below its derived upper limit. 

We discuss the origin of the [Ne n] and [N n] lines from 
diffuse low-ionisation gas in Sect. [5] and the possible origin of 
the [O m] emission from diffuse high-ionisation gas in Sect. [5] 

We also find that %// contributes very little to the [O i] emis- 
sion (3%), as expected since the atomic oxygen is a tracer of 
the neutral gas. ^hii shows very little contribution to the [C n] 
(<1%) and [Si n] (11%) lines as well. The energy required to 
create C ++ ions is 24.4 eV, close to that of S ++ , therefore when 
the radiation field is sufficiently hard, as is the case inside the 
H ii region, C is doubly rather than singly ionised. 

4.2. The dense photodissociation regions: ^pdr 

We subsequently investigate the properties of the neutral gas and 
compare model predictions to the [C n] and [O i] 63 and 145 fim 
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lines. To this end, we extend %// into the atomic phase and the 
outer envelope of the molecular phase, until the molecular frac- 
tion Fb/H reaches 99%. The modeled shell is fully covered by 
the gas. The ionised cloud, with average density of 10 2 8 cm -3 , 
transitions from the H n region into a dense PDR, which is where 
the [O i] originates. This transition between these two phases in 
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a constant pressure model translates into a jump of the gas den- 
sity at the ionisation front. With this model, the average density 
of the PDR is 10 5 1 cm" 3 . The radiation field striking the PDR 
is: Go = 3.5 x 10 3 Habing unit (1.6 x 10" 3 ergs cirT 2 s" 1 )- The 
heating of the gas in the PDR is dominated by grain photoelec- 
tric heating (85%), and there is a second order contribution from 
collisions with H°, He, and H2. The cooling is done by radiative 
de-excitation mainly of the [O 1] 63 /urn (65%), [Si 11] (20%), 
[C n] 157 jum and [O 1] 145 /urn (<7% each) lines. At such low 
A v <1, the [O 1] lines are optically thin. 

The [O 1] 63 and 145 fim absolute intensities are over- 
estimated by a factor of 10, while their ratio is reproduced within 
15% by this dense PDR component. Along with the [O 1] 63 //m 
line, [C 11] 157 /mi and [Si 11] 35 /urn are the brightest PDR lines. 
While the [Si 11] line intensity is over-predicted by a factor of 2 
by the dense model, the [C 11] line agrees with the observation 
within 10%. 

The consistent over-prediction of the two [O 1] lines may be 
explained by the fact that, in reality, the dense PDR does not 
cover the full dense H 11 region. Assuming that all of the [O 1] 
arises from a single dense phase (and not the diffuse gas, see 
Sect. 15.31 ). we reduce the covering factor o f the dense PDR to 
only 10%, which is a proxy for dense clumps dHunter et al.l 2001: 
lGrovesetalJ l2008). The PDR model with a covering factor of 
1/10 matches the observed [O 1] emission, but under-estimates 
the [C 11] emission by an order of magnitude. The main con- 
clusion from this is that the dense PDR component cannot ac- 
count for the bulk of the [C 11] emission observed in Haro 11, 
but only for 10%. This mismatch for the [C 11] is due to the rel- 
atively low critical density of [C 11] (2 x 10 3 cm" 3 ) compared to 
the high density reached in the PDR (npoR = 10 51 cm" 3 ), gener- 
ated by the pressure equilibrium and the drop of gas temperature 
between the ionised and neutral phases. Indeed, Cloudy PDR 
models with constant densities of ~10 35 cm" 3 and a covering 
factor of unity would reproduce the absolute intensity of the two 
[O 1] lines simultaneously, and yet under-predict the [C 11] by a 
factor of 2. Moreover, the H 1 mass of such models would exceed 
the upper limit of 10 8 M Q found bv lBergvall et all d2000t) . 

The small fraction of [C 11] that ori ginates from the dens e 
PDR is at odds with the finding of iBergvall et alJ ((2000), 
who estimate that more than 80% of the observed [C n] in 
Haro 1 1 comes from the neutral gas. They use stellar clusters 
of T eff ~ 35 000-40 000 K to account for the [O m] 88 fim 
emission, which contributes less than 20% to the [C 11] emis- 
sion, and they assume that the rest of the [C 11] (>80%) is associ- 
ated with PDR s . Then , they find npoR = 2 x 10 3 cm" 3 using the 
iKaufman et al.l d 19991) model, which assumes a constant density. 
The m uch higher density t hat we find is compatible with the find- 
ings of IVasta et alJ (2010) who analysed ISO observations with 
PDR models for several dwarf galaxies, including Haro 11, and 
find PDR densities between 10 4 and 9 x 10 4 cm" 3 , and Go be- 
tween 60 and 8 x 10 2 Habing. They also estimate that 10 to 60% 
of [C 11] comes from the ionised gas, which is in agreement with 
the diffuse ionised model contribution to the [C 11] line that we 
discuss in Sect. 15. 31 

We have also explored the influence of cloud geometry 
on the line predi ctions by using the PDR model Kosma-r 
(iRollig et alJl20 06). in which the illumination of the cloud is ex- 
ternal, simulating a spherical cloud structure. At corresponding 
values of n H (10 51 cm" 3 ) and Go (3.5 x 10 3 Habing) found with 
Cloudy, Kosma-r predicts ratios of [O 1] 145 fim/[C 11] ~1 and 
[O 1] 145/63 fim~Q.Q5. The [O 1] 145 //m/[C 11] model ratios are 
over-predicted, as found with Cloudy, compared to our observed 
value of 0.08. However, if we attribute only 10% of the [C 11] 



emission to the PDR, then [O 1] 145 yum/[C ii]pdr ~ 0.1 from 
Kosma-r, which is in line with our observations. The spherical 
geometry affects both [O 1] and [C 11] lines similarly to first or- 
der, and is not the primary reason for the discrepancy between 
observed and predicted [C 11] intensities. 

We prefer to impose the pressure equilibrium between the 
compact H 11 region and the dense PDR, because this gives phys- 
ical i nsight into the pr esence of dense but fragmented PDRs 
(e.g. lHunter et alJ l2001). Although the pressure equilibrium as- 
sumption is not valid on large scales (nor is the assumption of 
a constant density) and, in particular, is expected to break in 
molecular clouds, here we do not model the molecular phase. 
We are interested in the interface between the dense H 11 region 
and the PDR, the presence for which there is observational evi- 
dence. This model solution of the dense PDR with covering fac- 
tor of 10% is noted c ^pdr- Next, we will discuss possible ways to 
lower the density in the PDR to account for the [C 11] emission, 
by either combining a high density with a low density model 
(Sect. l53l l. representative of the multi-phase ISM of galaxies, or 
by including a magnetic field (Sect. lSTTT i. 

5. A model for the diffuse medium: ( / 

5. 1 . Need for a softer radiation field 

The major discrepancies between observations and model pre- 
dictions in the ionised gas are seen for the [Ne 11] 12.6 fim and 
[N 11] 122 fim lines, which are both under-predicted by a factor 
of ~5 in the dense H 11 region model ^hii, which is not unex- 
pected. With excitation potentials of 14.5 and 21.6 eV respec- 
tively, emission of N + and Ne + originates only from the ionised 
gas, from both the dense H 11 region and the diffuse medium. For 
[N 11] 122 fim, ^hii does not reproduce the observed value even 
when varying the 3 free parameters. Only r; n > 10 218 cm may 
agree (bottom left panel of Fig. [6j. The contribution of the com- 
pact H n region to this line is thus marginal; mostly because the 
radiation field produced by the young starburst is too hard and 
favors the ionisation of N ++ rather than N + . For the [Ne 11] line, a 
small change in the age of the burst would impact on its predic- 
tion, as discussed in Sect. 14.1.21 For the same values of density 
and inner radius, a model of age 3 Myr (instead of 3.7 Myr) 
would agree better with the observed intensity of [Ne n]. For a 
burst age of 3.7 Myr, higher densities (nH > 10 3 cm" 3 ) and in- 
ner radii (r; n > 10 21 6 cm) are required to reproduce the observed 
[Ne 11] (top right panel of Fig. |6|. However, such values of nn 
and ri n would no longer agree with the other ionic lines. In order 
to reconcile the modeled [Ne n] and [N n] with their observa- 
tions, a component with softer radiation field is required. This 
highlights the fact that one component is not enough the account 
for the emission of all the observed lines due to the presence 
of several phases in the ISM of galaxies with different proper- 
ties. Our physical picture of the ISM of dwarf galaxies is of a 
highly permeable medium, where some of the ionising photons 
from the central starburst escape the H n region and travel large 
distances. This means that in dwarf galaxies, larger effective vol- 
umes can be maintained in an ionised state, compared to those 
of dusty starbursts. 

To model this additional ionised component we have tried 
three different methods: (1) set the cloud further away (higher 
ri n ) under the same starburst conditions, and allowing the density 
to be a free parameter, (2) stop c €mi before it reaches the ioni- 
sation front and take the output spectrum as the input spectrum 
of a lower density medium, (3) opt for a softer radiation field by 
including a scaled local interstellar radiation field. However, all 
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of these methods fail to reproduce the observed line intensities. 
Either the radiation is still too hard (case (1) and (2) when stop- 
ping too early), or too soft (case (3) and (2) when stopping too 
close to the ionisation front). 

The radiation field that works best to explain the observed 
[Ne n] and [N n] lines is a representative stellar SED from 
the Kurucz library of relatively low temperature (see next 
subsection). This is supported by the fact that more evolved 
stars producing low energy photons are present in Haro 1 1 
(Mic heva et al.l 120101) but not accounted for in the compact 
model 'it? C&hii+^pdr)- This older population affects mainly the 
emission of [N n], [Ne n], and to a lesser extent [S in] . 



5.2. Properties of the model & 

We ran constant density models of diffuse ionised gas, stopped at 
the ionisation front, with stellar temperature between 20 000 and 
40 000 K. The best fit to the [N n] 122 /jm and [Ne n] emission 
is obtained for a stellar temperature of 35 000 K. Figure|9]shows 
how the [N n], [Ne n], [S m], and [C n] lines vary with density 
for several values of r; n . We refer the reader to the Sections ^. 1.11 
and 14.1.31 for an interpretation of the behavior of the lines. The 
range of possible solutions for the [N n] 122 //m emission are 
densities below 10 3 cirT 3 and rj n higher than 10 210 cm. [Ne n] 
also agrees well with these conditions. At T e ff of 35 000 K, 
[S m] also emits, but remains below its observed value. Very 
little [N m], [Ne m], [O in], and [S iv] are produced in this dif- 
fuse component, and no [O iv]. Although [N n] 122 /zm has a 
relatively low critical density, the density of the N + emitting 
medium is usually well constrained by the [N n] 205 //m line. 
Unfortunately, the [N n] 205 fim line is difficult to observe with 
PACS because of spectral leakage (see Sect. 12.21 ). The upper 
limit on the [N n] 205 yum line sets an upper limit on the ratio 
[N n] 205/122 fim of 0.75 which corresponds to a lower limit 
on the density of t he "diffuse" gas co mponent of 10 cuT 3 (e.g. 
iRubin et all [19941 lOberst et alJl2006h . The [N n] 122 /im line 
critical density sets an upper limit of 10 3 cm -3 . Models with 
densities between 10 and 10 3 cirT 3 work equally well for the 
"diffuse" medium in acounting for the [N n] 122 //m and [Ne n] 
emission. 



5.3. [C n] contribution from the model 9) 

Since c &pdr predicts very little [C n] 157 /mi, we expect a sig- 
nificant contribution of the diffuse medium to the C + emis- 
sion. C + is u sually found in the surface layers of far- UV illumi- 
nated PDRs dStacev etal .1119911: iNegishi et aT1l2001l) , but it can 
also come from the ionised gas, with a contribution up to 50% 
( Mad den et al.lll993t ICarral et al.l[T994t lHeileslll994t lAbel et alJ 
2005). With ionisation potentials lower than that of hydrogen, 
both C + and Si + can be produced by low energy photons or ex- 
cited by collisions with electrons or hydrogen atoms, depending 
on the degree of ionisation from the emitting medium. 

Figure [9] shows the predicted [C n] line intensities from 
the diffuse medium with a stellar temperature of 35 000 K, 
as a function of gas density. The contribution of the diffuse 
phase (ionised+neutral) is in filled circles, and the contribution 
of the ionised phase only shown by open circles. Calculations 
are stopped when molecules start to form and the fraction of 
H2/H reaches 10%. At low nn, the empty and filled circles 
are very close: a large fraction of the [C 11] can arise from 
the ionised phase, until [C 11] is collisionally de-excited by e~ 
(n c „, = 50 cm -3 ), and emits mostly in the neutral phase, its in- 
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Figure 9. Predictions for the model @ of the [Ne n], [C n], [N n], 
and [S in] line intensities (l m „dei), normalised by their observed 
value (l bs) as a function of density ne- Models are for a diffuse 
ionised/neutral gas, with constant density, stellar temperature of 
35 000 K, and are stopped when molecules start to form. The 
[S m], [Ne 11] and [N 11] lines emit in the ionised phase. The con- 
tribution of the ionised phase to the [C n] line is indicated with 
open circles. The vertical dashed lines indicate density values of 
the diffuse mo dels @i (n H = 10 cirT 3 ) and @ h (n H = 10 
defined in Sect. 



3 cm 3 ) 



tensity then increasing with density. At ri n = 10 22 2 cm, the [C 11] 
predictions decrease for ne > 10 2 cirT 3 , because the atomic 
phase is very thin as the material is cold and enters quickly into 
the molecular phase where the models stop. We disregard these 
models which do not predict enough [C 11]. The range of possi- 
ble solutions for the [C 11] 157 fim emission are ri n higher than 
10 2L4 cm, and are degenerate in ne- We can set a low-density 
case where ne = 10 cirT 3 - for example. Then ri n = 10 218 cm 
is the best fit, and [C 11] comes mostly from the ionised phase, 
accounting for ~45% of the observed value. We also consider a 



high-density case where ne 



10 3 cm" 3 , 



and ii, 



10 21S cm. In 



this case, [C 11] is mostly emitted in the neutral phase, account- 
ing for ~40% of the observed value. Both models fill a larger 
volume than the compact model 
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5.3.1 . Constraint from the H 1 21 -cm line 



The H i 21 -cm mass upper limit of 10 8 M from Bergval fet alJ 
(2000) is an important gas diagnostic of the diffuse atomic gas 
and brings a strong constraint that we can use to differentiate 
between the low-density and high-density case models. The H i 
mass from ^pdr is 10 7 M . While the models from Fig. [9] were 
stopped at a fraction of H2/H reaching 10% to compute entirely 
the ionised and atomic phases, most of those models exceed the 
H 1 mass upper limit. By stopping those models when the upper 
limit on the H 1 mass is reached, we find that the low-density 
case model contribution to the [C 11] line is unchanged since it 
comes from the ionised phase, while the high-density case model 
contributes to only ~20% of the [C 11] line intensity, since the 
atomic phase is stopped at lower Ay than previously computed. 

For the rest of this study, we refer to the low-density case 
model as 3%, and to the high-density case model as f^,. We can- 
not exclude the presence of a diffuse neutral gas, but in all cases, 
there is a prominent contribution to the [C 11] from the ionised 
gas. We estimate that ~40% originates in 2>i\ the rest coming 
from c {opbr and little from S>h- Moreover, the diffuse models do 
not contribute more than 15% to the [O 1] lines, which are reli- 
able tracers of the PDR. 



6. Influence from other possible components 

6.1. Magnetic fields 

Magnetic fields impact the model solution, dominating the pres- 
sure deep into the clo ud, and they are ex p ected to be impo r tant in 
star-forming regions dShaw et alJl2009h . lRobishaw et al.l {2008) 
found magnetic fields on the order of a few mG in ULIRGs, from 
OH Zeeman splitting. In local dwarf irregular galaxies, magnetic 
fields are found to be weak, of a few fiG to 50 nG in more 
extreme cases (e.g. NGC 1569, NGC 4214: IChvzv et al.ll201 it 
Keplev et al. 2 1 Ol l20TTh . No study on magnetic fields has been 
conducted for Haro 11, so we did not include them in our grid 
calculation. Nevertheless, in view of the fact that we have found 
that the density of the dense PDR plays a determining role in 
the predicted line fluxes, we investigate the influence of adding 
a magnetic field B mag of strength 1 fiG to 3 mG to the solution 
model ^ described in Sect|4] 

The magnetic field {B mag ) is expressed as a pressure term 
(Pmag) through the cloud by the following equation: 



mag 



B 

and B mag (r) = B, 

07T 



mag 



( 0)|^)\where« = 2/3(1) 



where r is the depth of the cloud, B mag (0) and «h(0) are the 
initial values of the magnetic field and hydrogen densit y. This 
density-magnetic field relationship from Crutcher ( 1999) corre- 
sponds to the case where gravity dominates over magnetic sup- 
port in the cloud collapse, and the magnetic flux is conserved. 
Figure [10] shows the predicted [C n] and [O 1] intensities with 
For low values of B mag , thermal pressure, ex- 



K «h Tg as , dominates and the densities in the 



increasing B mag . 
pressed as P gas 

PDR are high, with strong emission of [O 1], as seen previously. 
With increasing B mag , magnetic pressure starts to dominate over 
gas pressure and the transition into the PDR is smooth, resulting 
in lower densities in the PDR, and less emission of [O i]. The 
[C n]/[0 1] ratio increases with increasing B mag . The [C 11] pre- 
diction falls between a factor of 0.5 and 1 of its observed value 
for all models. Values of B mag up to 10~ 3 75 G predict enough 
[O 1] emission. For B mag > 10~ 3 ' 75 G, the models under-predict 



the intensity of the PDR lines, which would require extra input 
power, although we are still within the errors (factor of 2). 

However, when lowering the density, the size of the PDR 
layer increases, and the H 1 mass as well. With no magnetic 
field, we find M(H 1)= 10 7 M with a covering factor of -10%. 
For B mag > IQ- 4 G, we find M(H i)> 2 x 10 8 M for a cover- 
ing factor of unity. By increasing the magnetic field strength, 
we go from a picture of small dense clumps to a more dif- 
fuse extended medium, which again does not agree with the ob- 
servational upper limit on the H 1 mass. Moreover, high field 
strengths (> 100 fiG) may be found in the cores of star forma- 
tion, but are unlikely to hold on galaxy-wide scales. Models with 
Bmag > 10 1 4 G are therefore discarded and appear in the shaded 
grey area on Fig. |T0j The model with B mag of 10~ 4 G is noted 

Then, if we consider that all of the [O 1] 145 yum emission 
originates in the PDR, we can consequently scale the [C 11] in- 
tensity and we estimate the contribution of the PDR to the [C 11] 
line to be from 10% up to at most 50% in strong magnetic field 
cases. 
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Figure 10. Influence of the magnetic field on the PDR lines. 
Top: Model predictions of the [C n] and [O i] intensities (l mo dei) 
normalised by their observed value (lobs) as a function of mag- 
netic strength B„ mg . The effective density in the PDR is indicated 
for each B mag value. Bottom: Fraction of [C n] coming from the 
PDR, scaled to the [O i] 145 pan prediction, assuming that all of 
it originates in the PDR. Models in the shaded area are not valid 
with respect to the upper limit on the H i column density. 



6.2. Origin of the [O m] 88 pm line 

Although emission from the [O m] 88 pm line is expected in 
starburst galaxies, recent Herschel obse rvations show that it i s 
exceptionally bright in dwarf galaxies ( Madd en et al.l [2012b). 
In Haro 11, it is the brightest of all MIR and FIR lines. 
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With an excitation potential of 35.1 eV and critical density, 
n cr u= 5 x 10 2 cirT 3 , this line is expected to originate from a high 
ionisation and relatively low density medium. Effective temper- 
atures, T e ff > 40 000 K, are usually needed to expl ain the bright 
TO m ] emission in star-forming irregular galaxies (Hunte r et alj 
2001). When considering both models c €hu and @, we find a 
good match between model predictions and observations for all 
ionic lines considered, although the [O ra] 88 /mi line is under- 
predicted by the models by a factor of ~2. ^hii accounts for 60% 
of the [O m] 88 /mi emission, and & does not predict any [O m] 
emission; the ionising photons are not energetic enough. This 
factor of 2 is within the uncertainties of this study, however we 
discuss the possible origin of the [O in] 88 /mi line from a highly 
ionised more diffuse medium. Although the [O in] emission is 
explained by ^hii in the extreme case of Haro 11, this may not 
be the case in other metal-poor dwarf galaxies. 

In the star-forming region N ll-B of the Large Magellanic 
Cloud. iLebouteiller et all (1201 2l) find bright [O hi] 88 /mi emis- 
sion that originates from spatially extended high-excitation dif- 
fuse ionised gas, and can be modeled by O stars distributed 
across the region. In the case of Haro 11, the density in ^hii is 
too high to fully reproduce the observed [O ra] 88 /mi emission. 
As shown in Fig. [6] lowering the density in %// to densities 
< 10 2 cirT 3 would better reproduce the observed [O in] 88 /mi 
intensity. 

Given the two components that we have modeled ^hii and 
3l, we may consider a picture in which c €nu does not have a 
covering factor of unity, but is actually porous. This would allow 
some fraction of the ionising photons to escape the compact H n 
region and travel further away in a lower density medium. Such 
a configuration could explain the [O ra] 88 /mi emission, but 
would also lead to over-predicting several lines. In particular, in 
our case, predicted intensities of the [S iv], [O iv], and [N ra] 
lines would be a factor of ~2 higher than their observed values. 



6.3. Origin of the [O iv] line and X-rays 

With an ionisation potential of 54.9 eV, [O iv] can only be 
ex cited by high ene r gy sources. Its origin has b een discussed 
in iLutz etal . (1998); Schaer er & Stasinskal (Il999h . It can come 
from AGN activity (X-rays), very hot sources or shocks. In this 
Section, we explore the role of X-ray photons and/or shocks. In 
dwarf starbursting galaxies, a hot and young stellar population is 
usually responsible for the [O iv] emission. This has been con- 
firmed by I SO observations of the dwa rf galaxies IIZw40 and 
NGC5253 (ISchaerer & Stasinskal fl999l Signatures of Wolf- 
Rayet stars in the optical and P Cygni profiles in the far-UV 
(O vi lines), revealing the presence of O su pergiants and a very 
youn g burst in Haro 1 1, were reported bv iBergvall et al.1 12006) 
and iGrimes et al.l (l2007i) . In the case of Haro 1 1, we find that 
the young starburst containing Wolf-Rayet stars can indeed ac- 
count for the [O iv] 25.9 /mi line intensity. The luminos ity of 
the [O iv] line is -0.006% of L T ir. ISpinoglio etlril d2012h have 
compared the IR/submm emission lines to the IR luminosity in 
local samples of AGNs and starburst galaxies. They showed that 
the [O iv] line is weaker in starburst galaxies, about one order of 
magnitude less intense than in AGNs, and that starbursts follow 
the relation L[o iv] ~ 10 Ljir in the high IR luminosity range, 
which is consistent with that observed in Haro 1 1 . The [Ne v] 
(97.1 eV) at 14.3 /mi is not detected in the Spitzer/IRS spec- 
trum. The upper limit in Table [2] is consistent with our model 
prediction. Its absence confirms the fact that there is no AGN 
activity nor intense hard X-ray emission in Haro 1 1 . The X-ray 
component discussed in Sect. l3.3.T1 in our model has little effect 



on the intensities of the ionic lines. Around the parameters of 
adding X-rays to the radiation field increases the intensity 
of the [O iv] 25.9 /mi line by only -1%. In the PDR, the X-rays 
have a moderate effect. They contribute up to 5% to the [C n] 
intensity, and -10% to the [O i] intensities. X-rays play a second 
order role and the physics of the neutral gas is dominated by the 
FUV photons (PDR), not the X-rays (XDR). 

Shoc ks are also not our prefered explanation for the [O iv] 
emission. lO'Halloran et al.l ((2008) find that the emission of the 
[Fe n] 26.0 /mi line in low metallicity BCDs may be shock- 
derived or can result from a larger abundance of iron in the gas- 
phase as it is less depleted onto dust grains at low metallicity. 
However, the [Fe n] 26.0 /mi line is barely detected in Haro 11. 
The observed ratio of [Fe ra]/[Fe n] is ~0.005, which also tends 
to rule out shocks as a major heating process here. 

7. Multi-phase build-up of Haro 1 1 

We summarize the contributions of the different ISM phases that, 
when put together, account for the global MIR and FIR fine- 
structure line emission observed in Haro 1 1 (Figure fTTT l. The 
panel on the left is an illustration of this multi -phase build-up: the 
central starburst is surrounded by our main model c € composed 
of a compact H n region (red), with adjacent dense PDRs 
&PDR (blue). The volume around this compact nucleus is filled 
by diffuse ionised/neutral gas & (yellow). We also add a low fill- 
ing factor component of warm dust close to the starburst (black 
dots), the need for which we discuss (Sect. l7Tl i. This scheme is 
the global picture (not to scale) that results from this modeling 
study. We remind the reader that the models ^ and & were com- 
puted separately. The panel on the right is a histogram indicating 
the contribution of each model to the line intensities, with iden- 
tical color-coding. All 17 lines considered are reproduced by the 
3 models ^pdr, and & within a factor of 2. The model pa- 
rameters of these phases and line predictions are also listed in 
Table [5] In view of the results from Sect. 15.31 and 16.11 in which 
the upper limit on the H i mass limits the contribution from the 
diffuse neutral medium and from magnetic fields, the models &h 
and are not taken into account for our final build-up of the 
ISM of Haro 11. The models we combine are: 

- a compact H n region that dominates the emission of the ionic 
lines "tomi, 

- a dense PDR of coverage ~10% ^pdr, 

- a diffuse low-ionisation medium 

We analyse results from this build-up in terms of dust emission, 
extinction, and mass budget. 



7.1. The global SED 

The proposed geometry also implies a prediction for the dust 
emission coming from the different ISM phases of Haro 1 1 . We 
compare the predicted SED of the multi-phase build-up with the 
observations of Haro 1 1 in Figure [12] The goal here is not to 
fit in detail the dust emission but rather to reproduce the global 
shape of the SED because it is sensitive to the temperature of 
the dust, which is one of the main predictions of the multi-phase 
model. The SED of the compact model H (both c €uu and ^pdr) 
peaks at ~40 /mi (Fig.[T2l. as in the observed SED. In particular, 
ffpDR is too dense to explain the FIR shape of the SED, which 
corroborates the presence of a more diffuse component. Indeed, 
the diffuse component matches the temperature of Haro 1 1 
in the FIR, although the longest wavelengths, beyond 250 /mi, 
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Table 5. Model parameters and line predictions. 



Quantity 


Model 




+ VPDR 


% 


@k 


Warm dust 


Total {a) 


Energy source 


Starburst 3.7 Myr 


Star 35 000 K 


Star 35 000 K 


Starburst 3.7 Myr 




n H [cm- 3 ] 


io 2 - 8 10 5 - 1 ( 10 34 ) 


10 10 


10 30 


io 1 - 5 




Tin [cm] 


10 214 10 214 


10 21 8 


10 218 


IO 198 




I([S iv] 10.5)/U S w 


0.84 (0.34) - 


0.06 


_ 


_ 


0.90 


I([Nen] 12.8)/U 


0.20(0.12) - 


0.53 


0.89 


_ 


0.73 


I([Nem] 15.6VU 


1.13(0.72) - 






_ 


1.13 


I([Sm] 18.7VU 


1.03 (0.59) - 


0.26 


0.16 


_ 


1.29 


I([S in] 33.5)/W 


1.04 (0.75) - 


0.61 


0.14 


_ 


1.65 


I([Si ii] 34.8)/U 


0.11 (0.07) 0.24(0.37) 


0.04 


0.24 


_ 


0.39 


I(Hu a 12.4VU 


1.28 (0.72) - 


0.37 


0.63 


_ 


1.65 


I([Fem] 22.9VU, 


0.49 (0.32) - 


0.44 


0.50 


- 


0.93 


I([Orv]25.9)/U 


1.18(0.87) - 






0.10 


1.28 


T/r a„ __-| o nn\/T 

l(LAr in] 8.99)/l i,, 


1.49 (0.61) - 


0.08 


0.06 




1.57 


I([Arn]6.99)/U 


0.12(0.09) - 


0.75 


1.49 




0.87 


I([C n] 157)/W 


0.01 (0.01) 0.07 (0.38) 


0.40 


0.19 




0.48 


I([Om] 88)/U 


0.59 (0.52) - 


0.02 






0.61 


I([Oi] 63VU 


0.03 (0.02) 1.11 (1.02) 




0.08 




1.14 


I([Oi] 145J/U 


0.03 (0.02) 0.95 (0.98) 




0.07 




0.98 


I([Nm] 57)/U 


1.41 (1.17) - 


0.13 






1.44 


I([Nn] 122VU 


0.29 (0.28) - 


1.46 


0.46 




1.75 



(a) The total is taken as the sum of models ^ , S> t and warm dust, which is what we consider for the final build-up in Sect. [7] although and ^ B 
may also alter the line predictions. 

(b) Ratio of the predicted intensity of the model over the observed intensity (I j, s ). Ratio values below 1% are not indicated. 
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Figure 11. Left: Multi-phase build-up of the ISM of Haro 11, composed of a compact H n region (%//) in red, a dense fragmented 
PDR (&pdr) in blue, a diffuse medium of lower ionisation {&) in yellow, and warm dust in the inner region (represented as dots 
close to the starburst). Right: Diagram of the contribution of each individual ISM phase to the global emission of the MIR and FIR 
lines. Red color corresponds to blue to &pdr, with possible enhancement in the [C n] (and [Si n]) lines when considering 
magnetic fields (hashed blue, and yellow from the diffuse ionised medium also with possible contribution from a diffuse 
neutral phase (hashed yellow, 



are still under-predicted. Perhaps the most surprising result is 
that the multi -phase model fails completely to reproduce the UV- 
optical extinction and the MIR emission, despite the choice of a 
lower cut in the grain sizes (Sect. 13.3.31 ). We discuss these two 
points in Sect. l7.1.1l and Sect. 17. 1.21 



7.1.1. Energy balance 

In the UV-optical wavelength range, the total outward luminos- 
ity (Fig. |T2j bottom panel, black line) is 80% greater than the 
observed luminosity. Has the input energy in the UV-optical 
wavelength range been overestimated? We are confident about 
the shape of the input starburst spectrum since the spectral line 
ratios depend on the hardness of the radiation field and are well 
fitted by our models. Additional confidence draws from the right 
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amount of input energy as the level of the NIR photometry, unaf- 
fected by the little extinction present in Haro 11, is well matched. 

The SED reflects the global energy balance of Haro 1 1 and 
shows that a fraction of the total input energy injected in the UV- 
optical wavelength range is not absorbed. The models ^hii and 
S>t are responsible for the over-prediction in the UV-optical. We 
find the average extinction of the model to be Ay = 0.15 mag. 
The spectrum is dominated by emission from the ionised model 
components: ^hii, of which 90% of the radiation escapes since 
the PDR covers only 10%, and These models are com- 
puted separately and work as MIR-FIR emission line diagnos- 
tics, very little affected by dust attenuation, but we neglect their 
interaction. We can estimate the reddening needed to match the 
level of the model sp ectrum to the o bserve d data points. Using 
the extinction law of ICalzetti et alj d2000l) . with R v - 3.3 set 
by the model, we find Ay = 0.9 mag, hence Ag = 1.2 mag 
and E(B — V) = 0.3 mag. The reddened spectrum appears in 
dashed dark line on Fig. W2\( bottom). From UV contin uum mea- 
surements and Ha/Hfi ratio, Bergvall & Ostlin (2002) estimate 
Ab = 0.7 mag. The optical image (Fig. [TJ shows a dust lane 
passing in front of knot B, where the extinctio n is known to 
be hig her (I Adamo et al .1120101: iGuseva et al.ll2012l) . lAdamo et all 
(120101) find that E(B - V) is 0.4 towards knot B and 0.1 towards 
knot C. 

On the contrary, in the MIR (5 to 40 pm), the predicted lu- 
minosity is twice lower than the observed luminosity. This emis- 
sion in the MIR can be attributed to warm dust that we model in 
the following Section. 

This suggests an interpretation in terms of geometry of the 
galaxy, in which a large number of ionising photons can travel 
far, but in the end, the different media (compact and diffuse) are 
well mixed together, yielding a significant dust column along 
most lines of sight. In particular, the warm dust component that 
we add close to the central starburst and emits in the MIR is 
likely responsible for the absorption of the UV-optical part of 
the ^hii spectrum which is the primary contributor to the total 
spectrum in this range. 

7.1.2. Warm dust component 

To account for the observed warm dust that our models do not 
produce in the MIR range (FigQ~2l top panel), we have examined 
two possible solutions that do not change the intrinsic properties 
of the dust grains (which are little known): 

1) change the du st distribution for a postshock distribution 
dJones etal.l [T996). where grains are affected by shocks, with 
preferential destruction of large grains, resulting in an enhance- 
ment of small grains with respect to large grains. This grain shat- 
tering results in a distribution with steeper slope and lower cut- 
off sizes. This distributi on is successfu ll y applied to model IR 
SEDs of LIRGs in lDopitaet all J20TI . iGalliano et all l200l 
used this distribution, that may have been generated by super- 
nova shocks, to model the dust in the dwarf galaxy NGC 1569. 
However we have been unable to satisfactorily fit the SED by 
changing the grain properties. 

2) add a warm dust component close to the central starburst. 
In this way we add to our picture a porous layer of dust parti- 
cles close to the central starburst, wh ich is not inconceiva ble for 
compact H 11 regions. Models from iGroves et alJ ([2008) show 
that the observed warm dust emission arises from hot dust em- 
bedded within the H 11 region. This is also s upported by spatially 
resolved observations of nearby galaxies. ICalzetti et al.1 ([2005) 
find that the 8 pm and 24 pm emission originate within the H n 
regions. 



We adopt a warm dust component, at a distance of 
r; n = 10 19 8 cm, with density ne = 10 cm" 3 and covering fac- 
tor 15% (orange spectrum on Fig[T2b. These parameters are cho- 
sen to reproduce the level of the SED in the MIR. In these con- 
ditions grains can easily survive. This component has little or 
no effect on the discussed emission lines because of its com- 
pactness. The photons are predominantly absorbed by the dust, 
not the gas, and dust collisions are responsible for the cooling. 
This warm dust component ma y be regarded as an equivalent 
of an ultra-compact H 11 region dDopita et al.ll2.Q03h . or the hot 
spots described in Sieb enmorgen & Krugel (12007I) . This com- 
ponent takes away 15% of the input luminosity from the H 11 
region model. The amount of necessary energy and the hard- 
ness of the radiation field can be adjusted by modifying simul- 
taneously the input luminosity, covering factor, and the inner ra- 
dius in the code. Integrating the spectrum from 3 to 1100 pm 
gives Ltir = 1.7 x 10 11 L Q , very close to the infere d dust-model 
value of 1.4 x 10 11 L Q from lGalametz e t al. (2009); Remv et al.1 
(120121) . 

This warm dust model matches the continuum level in the 
NIR-MIR, but produces a silicate feature in emission at 9.7 pm, 
which is not seen in the Spitzer/IRS spectrum. This discrepancy 
is likely due to the fact that we have modeled the compact and 
diffuse phases of the ISM independently. Foreground dust from 
the other components may in turn extinguish the silicate seen 
in emission. Another factor at play is the adopted abundance of 
silicate grains. Quantitative surveys with Spitzer of AGB carbon 
and silicate dust production rates in the Magellanic Clouds find 
that carbon dust injection rates domina te over those of silicate 
dMatsuura et al.ll2009h iBover et al.ll2012l) . There may be a small 
effect from foreground cold dust from the molecular phase of 
Haro 1 1 that we do not model, which would have low filling fac- 
tor and high Ay. 

7.1.3. FIR-mm wavelength range 

The mismatch between the observed and predicted SED be- 
yond 250 pm (Fig. [T2l > may be attributed to colder dust in the 
molecular phase. With no molecular gas constraint, we cannot 
include a cold molecular phase in our modeling. However, as a 
first approximation on the required dust properties, we can con- 
tinue ffpoR to higher Ay and stop the model when it matches 
the submm temperature. This way, we estimate that the missing 
emission requires high Ay~40 mag and agrees with a grain tem- 
perature of ~25 K, which would fit the SPIRE data to 500 pm 
but still under-predict the LABOCA 870 pm data (see Fig. IA.ll 
in Appendix|A|. The latter is rather related to the submm excess 
introduced in Sect. 1 1.1 1 This excess is modeled by a cold dust 
component in Galametz et al.l d2009l) . 

7.2. Optical lines 

Several optical lines, including Ha, HyS, [O m] 5007 A, have been 
observ ed towards the central visual peak in Bergvall & Ostlin 
(120021) . Because of the small aperture of the observations, 
4" x 4", an absolute flux comparison cannot be made, e xcept 
for Ho- which was fully mapped in Ostlin et al. (1999) with 
the ESO 3.6 m telescope. More recent observations were also 
performed with the VLT/X- shooter, in the wavelength range 
3 000-24 000 A, bv lGuseva etall d2012l) . In particular, they ob- 
served, in a slit of size ~ 1" X 11", the Ha, HJ3, [O n] 3727 A, 
[O m] 4959 and 5007 A, [S n] 6716 and 6731 A, [S m] 9069 
and 9535 A lines, in both knots B and C of Haro 1 1 . From these 
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Figure 12. Top: Spectral energy distribution (SED) of the incident radiation field striking the cloud and outward radiation from 
each component modeled by Cloudy. The dark spect rum is the total ou t ward radiation r e sultin g fro m multi-phase build -up, and 
arriving to the observer. Photometry data from NED. iBendo et alJ d2012l) . [Galametz et al.l (2009), and Remv e t alJ d2012l) . and the 
IRS spectrum are overlaid in grey. Bottom: Spectral energy distribution of the total outward radiation field and the reddened outward 
radiation field in dashed dark line. 



observations, they determined the metallicity of Haro 1 1 and el- and comparison of the model predictions to the observations is 
ement abundances. Optical lines are affected by dust extinction, very dependent on the adopted geometry. With Cloudy, we can 
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compare the intrinsic emission predicted by our model to the 
extinction-corrected observations. 

The ratio of the optical lines listed above to H/3 are all well 
reproduced by our Cloudy models, within 40%, with contribu- 
tion from both the H n region and the diffuse medium (Table [6]). 
The only line which poses a problem is [S n], which is over- 
estimated by a factor of 3 at 6716 A, and 1.5 at 6731 A. The 
Clou dy ratio of the t wo [S n] lines is thus 0.9. The Cloudy mod- 
els of IGuseva et alJ (120121) also over-predict the [S n] lines by a 
factor of 2, but reproduce their observed ratio of 1.5. This ratio 
is sensitive to the electron density of the medium probed, which 
in their case is a density of 10 cirT 3 . We find that our diffuse 
ionised model with density 10 cirT 3 does predict a [S n] line ra- 
tio of 1.5. Concerning Ha, we compare its absolute flux and find 
that it emits from both the H n region and the diffuse medium, 
with a total pre dicted luminosity o f 7 X 10 s L , close to the ob- 
served value of lOstlinetai] (fT999), L Ha = 8 x 10 8 L . 



Table 6. Optical line predictions. 



Line 


I/KH/J) 




Observations*" 1 


+ 


% = 


Total Model* 6 ' 


[O ii] 3727 A 


2.3 


1.4 


3.5 


2.3 


H/3 4861 A 


1.0 


1.0 


1.0 


1.0 


[O m] 4959 A 


1.1 


1.7 


0.08 


1.5 


[O in] 5007 A 


3.4 


5.2 


0.03 


4.5 


Ha 6563 A 


2.9 


2.6 


2.0 


2.9 


[N ii] 6583 A 


0.53 


0.32 


0.79 


0.52 


[S ii] 6716 A 


0.22 


0.07 


0.03 


0.07 


[S ii] 6731 A 


0.12 


0.08 


0.02 


0.08 


[S m] 9069 A 


0.15 


0.17 


0.13 


0.19 


[S m] 9532 A 


0.41 


0.42 


0.32 


0.47 



(a) Ra tio of the observed line intensity over H/3 from Guseva et al. 
d2012l) . The fluxes are extinction-corrected and averaged over the two 
knots B and C. (b) Ratio of the intrinsic line intensity over H/3 predicted 
by Cloudy. 



not include the observed submm excess. When model ing the 
dust submm excess as cold dust, Gala metz et alJ (12009) find a 
total dust mass of 2 x 10 7 M , which would result in a higher 
D/G mass ratio (~l/50) compared to our current gas modeling. A 
significant fraction of ISM mass would therefore not accounted 
for by our models, part of which should sit in the molecular 
phase that we have not modeled. CO(1-0) is to date undetected 
(Berg vail et al.l 2000). The state of a reser voir of warm molecula r 
gas is addressed in a forthcoming paper (Cormier et al. 2012b). 
We would like to emphasize that the determination of the cold 
dust mass is very sensitive to the model and geometry adopted. 
Both the cool dust interpretation of the submm excess and the 
model accounting for the missing FIR emission (Sect l7.1.3b are, 
without further onstaints on the cold phase, ad hoc means to fit 
the dust SED. In particular, the cold dust interpretation of the 
submm excess has little physical support as it is unlikely to en- 
counter such very large columns of (cold shielded) dust in dwarf 
galaxies, and because the submm excess do es not seem to corre - 
late spatially with the densest ISM phases ( Galli ano" et al.l201lb . 



Table 7. Mass budget in solar mass M, 



Model 


M(H n) 


M(PDR)<"» 


M clust 




3.0 x 10' 




7.2 x 10 4 


^PDR 




1.2 x 10 8 


3.2 x 10 5 


% 


5.5 x 10 s 


1.2 x 10 8 


1.7 x 10 6 


Warm dust 






8.5 x 10 3 


Total 


5.8 x 10 8 


2.4 x 10 s 


2.1 x 10 6 


Literature 


~ 10 9 (b > 


2 x 10 s (c) 


~ 6 x 10 6 (d) 



(a) M(PDR) is the sum of the H i mass an d the H 2 mass in the 
PDR envelope untra ced by CO. (b) iBergvall & OstlirJ (2002). (c) 
iBergvall etaU ((2000), where they also report M(H i) < 10 8 M Q . (d) 
This i s the dust mass with out including the cold component of 10 K 
from lGalametz et alJ ([2009). Their total dust mass is 2 x 10 7 M G . 



7.3. Mass budget 

The gas mass of each ISM component can be calculated using 
the density profile in the cloud and integrating over the volume 
of each slice inside the cloud. We then multiply the gas mass 
by 1.36 to account for helium, and use the D/G mass ratio of 
the model (2.64 x 10~ 3 ) to derive the dust mass. The masses of 
each component are summarized in Table [7] We find a mass of 
ionised gas (^hii+^i) of 5.8 x 10 8 M , and a mass of neutral 
gas (^pdr+^i) of 2.4 x 10 8 M . This neutral gas mass accounts 
for the H i mass as well as the mass of IT formed in PDR 
envelopes before the CO formation. This H? layer not traced 
by CO is referred as the "dark gas" dWolfire et al.l l2010h . The 
PDR-H2 mass is ~ 10 8 M and comes from ^pdr, while the H 1 
mass comes essentially from the diffuse model 2#i with M(H 1) 
~ 10 8 M fixed by its upper limit. As expected, the ionised phase 
is dominant in H aro 1 1 . The ionised g as and PDR masses are 
similar to that of IBergvall et al.l (|2000). We derive a total dust 
mass of 2.1 x 1 6 M , comparable to the value 6 x 10 6 M of 
iGalametz et all (120091) . 

These dust and gas mass inventories do not include model- 
ing of the cold CO-cloud. The missing FIR emission discussed 
in Sect l7.1.3l results in a dust mass of ~ 2 x 10 6 M , which does 



8. Conclusion 

We have presented SpitzerflRS and Herschel/PACS spectro- 
scopic observations of the MIR-FIR fine-structure cooling lines 
in the starburst low-metallicity galaxy Haro 1 1 . Investigating the 
nature of the different gas phases in this galaxy, we have modeled 
the ISM phases of Haro 1 1 step by step with the spectral synthe- 
sis code Cloudy. The major model results can be summarized as 
follows: 

1. The IRS and PACS spectra show very bright MIR and FIR 
fine-structure cooling lines, the brightest of all being the 
[O 111] 88 //m line. The galaxy is undergoing young active 
star formation, as traced by the [O m] 88 //m line (35.1 eV), 
which is about 3 times brighter than the classical neutral 
gas tracers [C 11] 157 /mi and [O 1] 63 /mi. High line ra- 
tios of [Ne ra]/[Ne n], [S iv]/[S m], and [N m]/[N 11] also 
trace the hard radiation field. We find iij/Ltir = 0.1%, 
(L [C ii]+L [0 i])/L TIR = 0.2%, and altogether, the MIR and FIR 
lines represent 1.2% of Ltir. 

2. We model our exhaustive dataset of 17 fine-structure lines 
with Cloudy and find the need to describe Haro 1 1 with these 
component phases: 

- a dense H 11 region, illuminated by a young starburst, from 
which most of the ionic lines originate, 
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- adjacent dense PDR of low (~10%) covering factor, 

- an extended diffuse low-ionised/neutral medium, 

- a porous warm dust component close to the stars that ac- 
countes for the elevated MIR continuum. 

3. Emission from the [N n] 122 fim and [Ne n] lines cannot be 
reconciled with the compact H n region but are rather asso- 
ciated with a more diffuse low-ionisation medium. 

4. We find that 10% of the [C n] is associated with the dense 
PDR, and up to 50% when lowering the density. At least 50% 
of the [C n] arises from a more diffuse ionised medium. The 
low extinction and low metallicity of Haro 1 1 makes its ISM 
more leaky and the emitting region of [C n] larger. The [O i] 
emission is fully associated with the PDR. The observed 
[C n] luminosity is comparable to that of the [O i] 63 /im 
line. In the dense component we derive L[c ii]/L T i R ~ 0.05% 
and (Lp n]+L[o i])/L T i R ~l%, which is a common measure 
of the gas heating efficiency in the neutral gas, and in the 
diffuse component we find ii]/Ltir ~ 0.1%. 

5. The only line that we do not completely reconcile with our 
model is the [O in] 88 //m. Our compact H n region model 
accounts for 50% of the emission of [O in]. Our prefered 
explanation for the intense [O m] emission is the presence of 
yet another additional diffuse component, filling the volume 
around the porous H n region, and where the gas is heated by 
escaping UV photons. 

6. We estimate the mass from each modeled phase. We find 
a mass of ionised gas of 6 x 10 8 M , a PDR mass of 
2 x 10 8 M , and a dust mass (without sub mm constraints, 
because of the submm excess) of 3 x 10 6 M . The ionised 
gas mass is larger than the atomic gas mass. 

7. Finally, in terms of structure, the ISM of Haro 1 1 appears to 
be mostly filled with extended diffuse gas. Our simple pic- 
ture can reproduce the observations by a galaxy with a fill- 
ing factor of diffuse neutral and ionised gas of at least 90%, 
a dense H n region component of filling factor ~0.2% and a 
PDR component of filling factor <0.01%. A radiative trans- 
fer model which takes into account clumpy source and gas 
structures would provide a more realistic picture of Haro 1 1 . 
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Appendix A: SED including the cold molecular 
component 

We model the SED beyond 250 fim by stopping the dense PDR 
model ffpoR at larger Ay so that it predicts cold dust emission 
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Figure A.l. Same as Fig. Q~2](fo/? panel) except that the PDR component is continued to higher Ay in the molecular regime to 
the observed submm photometry data. The covering factor of this molecular component is of ~5%. 



